CHEMICAL ENGINEERING SCIENCE 
GENIE CHIMIQUE 


JUNE 1956 


Hold-up and pressure drop for water irrigating ‘‘ non-wettable’’ coke 


G. C. GarDNER 
Development Dept., Imperial Smelting Corporation Limited, Avonmouth. England 


(Received 20 June, 1955) 


Abstract—Simultaneous measurements of the liquid hold-up and of gas pressure drop, for 
beds of “ non-wettable ” coke (contact angle about 90°) irrigated by water counter-current 
to an uprising gas stream, were made for water rates from 0-8 to 12ft.*/hr. ft.” of column cross- 
section and for pressure drops up to 4° W.G./ft. height of column. Coke size ranges were 
and 

It is shown that the hold-up can be separated into two quantities —- that assumed to be 
in «a “ virtual state of rest ’’ and that flowing over the packing surface. The first part consists 
of a quantity //,, the “ static '’ hold-up, which will not drain from the packing when liquor 
and gas flow are stopped; a quantity //», which must be added to Ig when liquor flows, but 
which is independent of the liquor rate, and a quantity K, h" which indicates that some liquid 
is held in a state of rest by the pressure drop 4. The contribution to the hold-up from liquid 
flowing over the packing surface is expressed by means of a flow equation for the case of zero 
gas velocity, but the equation must be moditied to account for the retardative or dispersive 
effect of the gas on the liquid when the gas flows. The contribution of flowing liquid to the 
hold-up is shown by the last term in the complete hold-up equation for the experimental system. 


os 
H =Hgs Up K,h" + Ky Ve") 
The exponent r is constant for all packing sizes and equals 7-4. J/g, Ip, K,, Kg, Kg and n all 
vary with packing size but, with the exception of Hg, none can be correlated with that variable, 
as each shows a reversal of trend in passing from the smallest coke size to the largest. This 
anomaly may not be experienced for wetting systems. 

The pressure drop is correlated in terms of the true gas velocity. It is found that there 
are two gas velocity ranges, in the lower of which the liquor rate has no direct effect upon the 
pressure drop and in the upper one it has only a slight but definite effect. 


Lower velocity range =k, 
Ww 


w 


k, and ky are empirically correlated in terms of packing size. 

In conclusion, it is emphasized that the correlations in this work were made possible by 
using the true gas velocity, based on the free space in the packing, rather than the superficial 
velocity, and it is suggested that the mechanics of packed columns will never be fully understood 
until this variable finds greater use in theoretical investigations. 


Résumé— Les auteurs ont procédé & des mesures simultanées du hold-up liquide et de la perte 
de charge pour des lits de coke non mouillables (angle de contact d’environ 90°) arrosés d'eau 
circulant & contre-courant d’un gaz ascendant. Ces mesures ont été faites pour des vitesses 
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de l'eau allant de 0,8 a 12 ft. br. ft.? de section droite de la colonne et pour des pertes de charge 


allant jusqu’é 4” W.G. ft. de hauteur de colenne. 


Les dimensions du coke étaient de Poet Pal 


Les auteurs ont montré que le hold-up comprend deux quantités ; la premiére qui correspond 


Hg représente le hold-up ~ 


a un “ état virtuel de repos.” le deuxiéme a Pécoulement sur la surface du garnissage. 


statique.” la fraction liquide qui est retenue dans la colonne 


ne fenctionnant pas, I] p sajoute a 11, quand le liquide s'écoule et est indépendant de la vitesse 


du liquide, K, 
de charge h. 


” 


La part du hi ld-up, doe au liquide sécoulant sur la surface du gurnissage s‘exprime au 


indique qu'une partie du liquide est retenue en état stationnaire par la perte 


moyen d'une équaticn d'écoulement dans le cas d'une vitesse gazeuse nulle, mais cette équation 
peut étre modifi¢e pour tenir compte de leffet du gaz. 


Pour le systéme en expérience, la contribution de lécoulement liquide est mesurée par le 


dernier terme de léquation complete. 


= He + Up 4 


K, h,, 


L’exposant r est constant pour toutes les dimensions de garnissage et égal a 7,4. [/15. Hp. 


Ke + Kaba’ 


K,. Ky, Ko et n varient tous avec la dimensicn du garni , mais a l'exception de //,. on ne 

peut en relier aucun a cette variable car chacun s‘inverse de sens en passant dun coke de trés 
petite dimension a un coke de trés grande dimension, Cette anomalie ne peut pas étre étudiée 


a deux domaines pour cette vitesse : 


INTRODUCTION 


This work is part of a programme designed to 
give an insight into the manner of flow of slag 
over coke counter-current to an uprising gas 
streem, and later to indicate the heat and mass 
transfer rates to be expected between the slag 
and the gas. The informetion in this, the first 
part of the work, is solely concerned with the 
dynamics of flow of the two phases, and thus 
presents correlations of pressure drop and 
hold-up in terms of the gas and liquid flow rates. 

It would be extremely diflicult to carry out 
accurate experiments on the high temperature 
slag-coke system, and it was decided to use a low 
temperature model of water flowing over coke 
treated with a silicone fluid to make it non- 
wettable. In this way allowance was made for 
the principal difference between the slag-coke 
system and_ those investigated by 
other authors in connection with packed absorp- 
tion and distillation columns. <A difference, which 
may be regarded to be of secondary importance, 


systems 


expérimentalement pour des systemes mouillants, 


Les auteurs relient la perte de charge en termes de vitesse vraie du gaz. Ils trouvent qu’il y 
dans le plus faible, la vitesse de la liqueur n'a pas d'effet 
direct sur la perte de charge, dans le plus fort, leffet est faible mais détini. 
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was that the liquor flow rate range studied was 
0-8-12ft.*/ ft.*hr., whereas other investigators have 
been more interested in flow rates above 
10ft.” ft.2 hr, for in absorption and distillation 
columns the flow rate is easily controlled by 
recirculation or reflux, and high rates are normally 
desired for the best use of the packing surface. 


History 


Previous work on the fluid dynamics of 
packed columns can be divided broadly into two 
groups. 

1. Work making qualitative observations of 
the flow conditions, together with empirical 
correlations of pressure drop in the gas 
stream, hold-up of liquor and the transfer 
coefficients per unit volume. 


2. Work which has been interested in 
determining the actual liquor surface 
available for the transfer processes or the 
actual wetted surface of the packing. 


Of the first group the work on hold-up 
and pressure drop exemplified by Macn [7], 
Ucnipa and Fvgsita [13], and 
Wertss [1] and Jesser and Excoin [5] is of 
immediate interest. Qualitatively it was 
observed that at zero gas and liquor flow 
there was a static hold-up which would not 
drain away and which resided at the points 
of contact between the packing elements. 
With liquor flow, but no gas flow, the 
liquor flowed over the packing surface, 
and it was only at the higher rates that 
any great proportion of the liquor left 
the packing to pour through the gas, before 
recontacting a solid surface. The data 
for hold-up at zero gas flow has recently 
been correlated by Orake and Oxapa [10], 
who showed that the hold-up increases as 
the superficial liquor velocity to the power 
of 0-51 or 0.676, depending upon a modified 
Reynold’s Number. This correlation will be 
discussed in more detail later. 

With the liquor rate held constant, the hold-up 
showed no appreciable increase with respect to 
the superficial gas velocity until a point was 
reached in the loading range, or just before 
flooding, when it rose sharply. A similar effect 
was observed for the pressure drop, which 
increased steadily to approximately the same 
critical point, and then rose sharply to flooding. 
No attempt was made to correlate hold-up in 
terms of pressure drop. On the other hand at 
least two sets of investigators, Uchida and 
Fusita [18] and Reev and Fenske [11], have 
formed correlations for the opposite effect of 
hold-up on pressure drop for atmospheric and 
vacuum systems respectively, but their presenta- 
tions are such that no accurate reassessments 
of the data could be made to reverse the relative 
importance of the variables. Both sets of 
investigators employed the hold-up to account 
for a decrease in free volume of the packed 
column, but neither plotted their data directly 
in terms of the true gas velocity. Recently one 
other investigation, Lerner and Grove [16], 
has emphasised the importance of hold-up in 
modifying the gas velocity. A simple correlation 
was made for the flooding point, which implies 


Hold-up and pressure drop for water irrigating “ non-wettable " coke 


108 


that flooding is incipient when a critical true gas 
velocity is reached. As the hold-up is virtually 
constant for a given liquid rate for gas velocitics 
up to the flooding point, the theory was easily 
tested from published data. 

The second and more recent set of investigators 
has started from the viewpoint that, although 
the transfer coefficients per unit volume are the 
most useful to employ in practice, they should 
be separated into transfer coefficients per unit 
area and the transfer areas per unit volume. 
The problem is to measure one or other of the 
quantities, but unfortunately the various methods 
employed have given very diverse results. 
Mayo, Hunter and Nasa [8] measured the 
extent of staining of paper rings, GrimLey [2] 
measured the electrical conductivity of the 
liquid over a_ certain column height, VAN 
KreveLen, Horriszzer and van Hooren [14], 
Weisman and [15] Snutman and 
DeGourr [12] and Ismino and Orake [4] 
measured transfer rates in packed systems and 
compared them with rates obtained from surfaces 
of known area. 

There are a number of objections to those 
investigations concerning transfer processes, 
some of which could have been overcome by 
preliminary hold-up and pressure drop experi- 
ments. Most of them used transfer correlations 
employing the superficial gas velocity, and no 
allowance was made for changes in the velocity 
when the void space was decreased by hold-up. 
This may be a significant factor even if the 
static hold-up alone is considered. IsHino and 
Orake made the almost certainly invalid assump- 
tion that the liquid velocity was the same as if 
the whole void space was occupied by liquid. 
Secondly, if the mass or heat transfer coefficients 
could be used to measure surface area, the 
momentum transfer coeflicient could be used 
in the same way, by comparing pressure drops, 
and the error made by vaN KREVELEN in assum- 
ing that the total surface area exposed to the 
gas was equal to the surface area of the dry packing 
could at least be reduced. Lastly, interpretation 
of the hold-up results in terms of known flow 
equations might allow some conclusions to be 
reached as to the condition of the fluid in the 
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column, which would be useful in explaining the 
other transfer data. It is for these reasons that 
the experiments in this work have been performed 
before any mass-transfer experiments have been 
made. Although the system investigated may 
be regarded as a special one, it is believed that 
the results have implications of general value 
to irrigated packed columns, and may find 
applications outside the slag-coke system which 
it is primarily required to elucidate. 


THEORETICAL 

The true gas velocity V3, which is the gas 
volume per unit time per unit packed cross- 
section divided by the fractional free volume 
of the packing, taking into account hold-up, is 
used throughout this work in preference to the 
superficie! gas velocity, which is merely the gas 
volume per unit time per unit packed cross- 
section. 


1. Pressure Drop. 


The pressure drop through dry packed columns 
has been the subject of investigation by numerous 
authors. For the present purpose the correlation 
given by Mort [9] is suitable as it uses the 
surface area of the packing as a variable, thus 
allowing the calculation of the total exposed 
surface area in an irrigated column from pressure 
drop measurements. 

For single phase flow through a packed 
column the friction factor ¢, is plotted versus 
the Reynold’s number Re where ¢ and Re are 
defined by the equations 


(1) 
pol f/k)(1—e) 
or, if we desire to calculate the surface area, by 


Fig. 2 shows the various published data for 
spheres replotted after this pattern by Mort 
with the shape factor (f k) proposed by Hrywoop 
[3] equal to 6. 

When irrigated columns are considered, the 
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specific surface S becomes an unknown variable 
which can be correlated empirically in terms of 
the liquor flow-rate and the gas velocity. The 
preferable method of doing this is to correlate, 
first of all, the pressure drop empirically in terms 
of the same two variables and to substitute 
into an equation fitting the ~¢ — Re plot. The 
pressure drops are correlated in this work but 
the process has been taken no further, although 
the required information is presented, as the 
total surface area will only become of value when 
information on heat and mass transfer is obtained. 


2. Hold-up. 

An equation for the total hold-up is to be 
derived, which, though largely empirical in its 
formulation, will contain the postulation that 
the hold-up may be divided into two parts. 
One part is in a virtual state of rest and resides 
at places which have a natural retentive effect 
upon the liquid, such as points of contact 
between packing elements. The state of rest 
docs not mean that there is no liquid flow through 
this part, but that the flow is so slow that its 
volume is governed by the other effects such as 
the surface tension of the liquid. The second 
hold-up quantity is the flowing part, the volume 
of which will depend upon factors such as the 
amount of liquor flowing, its cross-sectional «r+ 
of flow, its perimeter of flow and the kinematic 
viscoscity. 

The hold-up in a virtual state of rest includes 
the static hold-up H.,, or that hold-up which will 
not drain from the packing when both gas and 
liquid flow are stopped. To H, must be added a 
small quantity H» due to flowing liquid passing 
in and out of the sites of static hold-up. H,» may 
have some dependence upon the amount of 
liquid flowing, but its quantity is so small that 
the analytical procedure cannot show any 
variation. When gas is flowing the pressure drop 
may increase the volume of sites containing 
liquid in a state of rest or may create new sites 
so that a further hold-up volume K,h”" must 
be added to Hg and Hy», where K, and n are 
empirically determined constants applicable to 
one system. Thus the total hold-up in a state of 
rest becomes 
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Hy + Hy + K, h" 


To obtain a logical form for the flowing part 
of the hold-up, consider liquor flowing in films 
or trickles over a vertical surface. The cross- 
sectional area of flow per unit packing cross- 
section is a, the perimeter of flow per unit packing 
cross-section is b and L is the liquor volume per 


per unit height 
my b\™ b L* 
L PL 


ax | (8) 
where the bracketed part of the proportionality 
is a suitably chosen Reynold’s number necessary 
in the friction factor. Now for a_ particular 
packed system it would be expected that the 
perimeter of flow would be solely dependent 
upon the area of flow, and it is assumed that the 
empirical relation 


(4) 


holds, where p would vary from zero, for the case 
where liquor was flowing over the whole packing, 
upwards. 

Rearranging [8] and substituting [4] 


b xa’ 


(5) 


m 2-m 
PL 


If the effect of both kinematic viscosity and of 
liquor flow rate the hold-up is known, values of 
m and p can be calculated to conform with this 
equation. This investigation dves not study 
variation of kincmatic visocosity, but reference 
may be made to the work of Orakre and OKxapa 
[10} who recently correlated all published data 
for wetting systems, to indicate how equation 
(5) may be applied, though not to give any 
quantitative information for the non-wetting 
system. For higher liquid flows but with no gas 
flow, where the effect of H, would be small, 
they found. 

Hy « 

BL BL 
where Hyp is the ‘dynamic’ hold-up or that 
which will drain from the packing. To make 
equations (5) and (6) consistent m = 0-467 
and p = 0-500, which implies that, although 


(6) 


the packing may appear completely wetted at 
high flow-rates, the perimeter over which the 
liquid is effectively flowing is still increasing 
since p is far from zero. 

Equation (5) does not allow for pressure drop 
in the gas phase. If this is included 


m 
+m) 


3—pi(l+m) 


There is a second and more important effect 
of gas flow, in that it retards the liquid by 
frictional effects or it disperses it into thinner, 
or more and smaller, trickles. This is accounted 
for empirically in the last term K, Vg, in the 
complete hold-up equation 


+ K, VG] (7) 


The first three terms in the right hand side 
represent liquid in a virtual state of rest, and the 
remaining term represents liquid in a flowing 
condition, 

APPARATUS 


The apparatus is illustrated in Fig. 1. The 
shaft, 24° tall and 7-5" x 8-5” in cross-section, 
was constructed from tinplate. Coke was 
supported at the base of the shaft by a grid 
of parallel 4," wires with § gaps between wires, 
and the top of the coke was enclosed by 4” 
wire netting to prevent the tendency of the coke 
to lift under the higher pressure-drops. Water 
was supplied from a head tank through a valve, 
which could be turned on or off rapidly, to a 
distributor containing either 33 No. 1 hypodermic 
needles or 33 3°’ lengths of 1 mm bore capillary 
tubing. Water was discharged about 8” above 
the coke surface. 

The base of the shaft was enclosed by a wind 
box. Air, previously metered by a sharp-edged 
orifice, was introduced into the base of the wind 
box and was distributed over the cross-section 
by a plate containing 64 3’’ D holes. Water 
draining from the coke bed was collected above 
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Water leaving the apparatus passed through 
HEADTANE a water seal, and thence to a drain or to a measur- 
ing cylinder. The gas pressure was measured 
just below the packed section by a sloping or a 
vertical manometer, according to whether the 
pressure was less or greater than 1" W.G. 


WATER 
DISTRIBUTOR EXPERIMENTAL PROCEDURE 


1. Preparation of the coke. 


A batch of Nantgarw metallurgical coke was crushed 
and graded to the size ranges 3° to 4° mesh, $” to 3” 
mesh, and 3’ to 1° mesh. Suitable volumes of each size 
range were dried in a steam oven then soaked in carbon 
tetrachloride and dried again. Finally the coke was 
soaked in a 5 volume °, mixture of Midland Silicone’s 
M.S. 200 silicone fluid and carbon tetrachloride, and the 
latter was evaporated. Water had a contact angle of 
about 90° with the resultant coke. 


2. Dumping and voidage of the coke column. 


a, Approximately an 18° depth of a weighed quantity 

TROUGHS WATER TO DRAIN of coke was dumped in four lots into the shaft of the 

OR CYLINDER apparatus, with vigorous banging between each lot to 

AIR OISTRIBUTOR F — — — ensure a stable coke packing which would be maintained 

a for a long period of time. The top of the coke was held 
FROM FAN down by the wire netting previously mentioned. 

Fic. 1 Diagram of apparetus. The apparent density of each coke size range was 


measured by weighing and by water displacement, 
the distributor plate , by two sets of parallel and thus the fractional voidage of the coke bed could be 
troughs, indicated in Fig. 1, which sloped towards  cajculated after an accurate measurement of the coke 
a common water off-take. depth had been made. 


ao Yain- COKE 
o Vzin- Yain COKE 
«x Yain-11n COKE 


oe ~< R.A. MOTT = BED OF SPHERES 


° 
Bo a 


Wr THE FRICTION FACTOR 
° 


200 300 400 300600 800 1000 
Re THE MODIFIED REYNOLOS NUMBER ——= 


Fic.2 Modified friction factor versus modified Reynold's number, for dry coke beds. 
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3. Measurement of the static hold-up. 


A dilute sodium chloride solution of known strength 
was prepared and 2 litres at a time were emptied on to 
the top of the coke in such a manner that flooding con- 
ditions were obtained. The solution collected from the 
base of the column was analysed for sodium chloride, 
and, when it showed correspondence with the inlet 
solution, the column was left to drain for about half an 
hour. The procedure was then repeated with distilled 
water until no sodium chloride could be detected in the 
draining water. Analysis of the wash water and its 
volume gave the necessary data for calculation of the 
static hold-up. 

The process was repeated for several heights of packing 
so that end effects at the top and bottom of the column 
could be eliminated. 


4. Measurements of the dynamic hold-up and pressure 
drop. 

Hold-ups for one coke size range were all made with 
the same dumping of coke. Water at the desired rate 
was fed on to the top of the packing and the air was 
turned on. It was assumed that a steady state had been 
reached when the pressure at the bottom of the column 
and the air orifice differential showed no variation with 
respect to time. Then the water rate was checked, by 
timing the filling of a 2 litre cylinder with the effluent, 
and recordings were made of the air flow rate and the 
column back pressure. The hold-up was measured by 
turning off the water supply, and later the air flow, and 
catching all the water that drained. Previous measure- 
ments had been made of the hold-up of the system for 
various water flows with no coke present. These relatively 
small values were subtracted from the drainage figures 
to give the dynamic hold-up of the system. 

Results were checked by staggering the measurements 
with respect to air volume. No measurements were made 
with u pressure drop of greater than 6° W.G., which 
was near the flooding point, to prevent any danger of 
the coke shifting and its packing characteristics altering. 

Before the coke was wetted a series of measurements 
of the back pressure of dry coke were made. After wetting, 
the back pressure was measured with no water flow. 
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L = 3-99 ft.3/ft.*hr 
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L = 2-39 ft.4/ft.*her 


Ve 
2-32 
3-66 
447 
5-07 


5-70 
6-35 
6-49 
6-63 
6-85 


Table 1. Experimental Results for }'’-4"’ Coke. 


Hp 
0037 
0038 
0049 

0056 5-64 
0064 
‘187 -0081 
0116 
0150 
0185 


| 0039 


“249 
‘278 
324 


RESULTS 
All the experimental results of true gas velocity, 
pressure drop and hold-up are given with respect 
to liquor rate and coke size in Tables 1, 2 and 8. 
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1. Pressure drop. 

(a) Pressure drop through dry coke. All the 
results are plotted as a ¥ — Re plot in Fig. 2 
using a shape factor (f/k) = 10 as suggested by 
Heywoop [3]. The trend of the curve is the 
same as that for spheres but is slightly displaced. 
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Table 2. Experimental Results for }''-}'' Coke. 


L = 0-80 


L= /ft2hr L = 3-99 ft.3/tt2hr 


| L = 7-08 ft.9/ft2hr = 11-06 ft.3/ft.*hr 
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(b) Pressure drop through wetted coke with 
no water flow. All the results are plotted in Fig. 3 
versus the true gas velocity. They are compared 
in the same figure with the results for dry coke. 
There is a considerable increase in back pressure 
on wetting the coke, indicating an increase in the 
specific surface S. 

(c) Pressure drop with water flow. The results 
were divided into two velocity ranges and empiri- 
cally correlated, as shown in Figs. 4 and 5, 
by the equations 


Lower velocity range h (%) =k, V;'* 
Ww. 


Upper velocity range h (° = k, 
pw 


where the V, is the air velocity in ft./see at 
15°C, 

The constants are given in Table 4. It will be 
noted that, for the range of liquor flows inves- 
tigated, the 4’° — $’ coke had no upper velocity 
range. 

k, can be empirically correlated by the equation 
106 x 10+ 
k, with less justification can be correlated by 
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| 
| 
Ve | | Mp | Vo | & | Mp | Ve| h | My 
| | | | 
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886 848 | 7-20 | -227 | | 7-59 | -259 | 0480 
O14 | 8-75 | -315 = 7-36 | -216 | | 7-77 | -281 | 0531 
925 | -313 | | | 7-50 | 251 | mmm | 7-96 | 289 | 
7-77 «-262 | | 8-12 | | 0618 
| | | 822 289 | 0476 
| 
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Table 3. Experimental Results for }''-1" Coke. 


L = 0-80 ft.3/ft*hr | L = 1-60 | = 2-30 | L = 7-08 ft.9/f. "he |L = 11-96 ft.3/ft.*hr 


4.85 x 10° 
d 


k, = 


3. Hold-up. 
(a) Static hold-up. Table 5 shows the results 
for the three coke ranges. Empirically it may be 


written 

2 x 10° 
Hs = (8) 
but it is doubtful if this expression can be 2-1” 
extrapolated very far. 


| | 
Ve h Hp | Va | h Hy | Va | | Hy | Va | Hp | Vg | | Hp 
207 | 010 0032 | 1-99 | 10 | 0045 | 2-06 | 010 0048 | 2-46 | 014 | -0109 | 8-48 | | 0442 
2-51 | O14 | | 313 | | | 2-52 | O16 026 | | 8-67 | -284 | -0516 
2-92 0032 | 3-82 | 031 | 0048 | 0054 3-81 | 030 | O11] | 8-75 | | 0546 
351 | 026 426 | 039 | 0051 | 3-97 036 | 0055 | 485 | 053 | 0114 | 9-03 | 332 | 
OL. 4:36 0032 | 5-70 | 068 | 0054 | 476 051 530 | 063 | 0118 
956 472 048 | 0035 | 6-80 | -107 | 0072 | 5-47 | 068 | 0065 | 5-76 | -082 | 0149 
| 
5-60 069 0025 ‘131 | 088 629 | 092 | 0069 | 6-36 | -103 | 0172 
6-28 | 092 | 0048 ‘172 7-00 | -126 | 0085 | 6-85 | -120 | 0172 
685 | 099 | 167 0119 | 7-66 | -161 | -O118 | 7-95 | -187 | 0235 
7-60 | -145 | | 8-72 | 213 | -o169 | 8-28 | | 0155 | 8-85 | 221 | 0917 
8-08 | 165 | 0096 259 | 0228 | 8-85 | | 0198 | 8-91 | -261 | 0382 
8466 | 184 vo | 975 | -283 | 0264 | 9-62 | -279 | 0270 | 9-02 | -267 | 0394 
9-31 | 234 9-45 | | 0567 | | 
9-57 | -263 | 0220 “330 | 
| | | 9-68 | 348 | 0663 | 
0-00832 ons 
0-00430 0-000927 
0-00274 0-000666 
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Table 5. 


kK, Ky 


000031 : x 109 5.385 x 107° 
0.00271 284 x 10% 72400 x 


0.00176 1715 x 10% 2-347 x 


x ORY COKE (b) Dynamic hold-up with no air flow. No 
© DRAINED COKE direct measurements of this quantity were 
made and data had to be obtained by extrapola- 
tion. This was not difficult as the hold-up 


WATER RATE 


0-78 cuft/sqfthr 


ft. H,0/ft 


2 
3 
7 


. 
. 
. 


in COKE 


Yain-Y/zin 4 


= 
3 
alé 
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Af PL 
( be) PRESSURE DROP PER UNIT COKE DEPTH ——> 


1A. 
[*¥ain-1in Cone 

Vg TRUE GAS VELOCITY——> ft/sec aT 15°C 2 3 4 $6789 


Fico. 8 Comparison between pressure drop over dry and Vq TRUE GAS VELOCITY ——» ft/sec aT 15°C 
over drained coke. Fico. 4 Pressure drop for lower velocity range. 
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Coke size | | 
Range €p Is Hp K, | 
0-417 0-0321 
0-462 0-0163 
at 
|, 
— 
|| 
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curves with respect to gas flow are very flat 
for small gas flows, but the data may be slightly 
less accurate than the measured hold-up results. 

Equation (7) indicates that the dynamic 
hold-up (H — Hs) = Hp should plot as a straight 
line versus a power of the liquid flow rate. 
Fig. 6 shows that an equation 


Hp = Hy + 


gives an excellent fit. Values of H, and K, 
have been calculated from these lines and are 
given in Table 5. A reversal in their values 
passing from through to }’-1" 
coke will be noted. It will be the subject of 
discussion later. 

(c) Dynamic hold-up at constant air velocity. 
Referring to equation (7) the exponent of L 
is set equal to 0-8 and that for (1-h) arbitrarily 
equal to 4. The value of the latter exponent is not 
critical, and 4 is of the correct order of magnitude 


WATER RATE 
cuft/sqft he 


4 6 0 
Ve TRUE GAS VELOCITY ft/sec aT 1s °C 


Fic. 5 Pressure drop for upper velocity range. 


to fit the theory. The exponent n and K, in 
K, h”" were obtained by plotting the rearranged 
form of equation (7). 


(Hp — Hg) (1 — K,h"(1 
=-i,+— 


as shown in Fig. 7 for various constant gas 
velocities where K, is a constant for a constant 
gas velocity. Choosing n by trial and error gave 
parallel lines which are straight in accordance 
with the reasoning of this paper. Values of n 
and K, are given in Table 5. 

(d) The effect of gas velocity on the dynamic 
hold-up. Equation (7) was written 


H — Hy, —Hy — K,h") (1 
( Ss 1 ) ( K,| =K,Vg 


The left hand side of this equation was plotted 
versus VW, and K, and r were determined with 
surprising result that, although K, varied with 
coke size, r was constant at 7-4. This conclusion 
is illustrated in Fig. 8, which also shows the 
accuracy of the correlation for, although a 
residue function of the hold-up has been plotted, 
an excellent curve is obtained. A few experimental 
points lying well below the general trend are due 
to results obtained with low water flows where 
the error is naturally greatest. 


T T 
Yein-Yain 


° 


Hp OVNAMIC HOLD-UP AT ZERO CAS FLOW —= 
° 
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Fic. 6 Dynamic hold-up at zero gas flow in relation to 
water rate. 


lll 


Seer.e™. 
04025270 


G. C. GarRDNER 


Ya in-¥21n COKE 


T T 
Ve in- Vain COKE 


T T 
Vain -11n COKE 


= 
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Fic. 7 The dynamic hold-up in relation to the pressure drop. 


The final hold-up correlation for water flowing 
over non-wettable coke now reduces to 
H =H, + Hy + K,h" + 
The first three terms on the right hand side 
are, according to the present theory, hold-up 
contributions present in a virtual state of rest 
condition in the coke column. The final term is 
hold-up present in a state of flow on the coke 
surface. 


Vain ~/2in CORE 


© Yzin-Yain Cone 
© Vain-1in COKE 


° 


ke pee 


Ky 


-$ 
° 2 a 6 8 10 


Ve TRUE GAS VELOCITY ft/sec at is °C 
Fio. 8 Correlation of hold-up residuc. 


Hs can be expressed in terms of the coke size 
by the empirical expression 

2 x 10° 
— 

The variation of all the variables in equation 
(9) with respect to the coke size are given in 
Table 5. 


Hs 


Discussion 
1. Pressure drop and specific surface. 


The y — Re plot for dry coke allows the calcula- 
tion of the specific surface of the irrigated coke 
system if Heywood’s shape factor value (fk) = 10 
is assumed correct, or if we use a slightly different 
(f/k) value to make the coke results cover-plot 
the results for spheres. A_ better alternative 
would be to compare the specific surface of 
irrigated coke with that for dry coke, by using 
the ratio of these two values, when no assumption 
as to the shape factor need be made. All the 
required data for the estimation are in this 
paper. As a guide, the ratio of the specific 
surface of irrigated coke to that for dry coke 
appears to be in the range from 1-0 to 1-2 for 
i’-4" coke, from 1-1 to 1-6 for }'-}" coke, 
and from 1-2 to 2-0 for }’’-1" coke. 

The pressure drop is correlated empirically 
in terms of liquor rate, true gas velocity and coke 
size. The use of the true gas velocity is important 
as it indicates a lower gas velocity range in which 
the pressure drop is independent of the liquor 
rate, and that even in the upper velocity range 
the effect of liquor flow, though definite, is only 
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slight. For the 47-4" coke there is no upper 
velocity range, but there is no reason why it 
should not appear at higher liquid rates, for it 
will be noticed that as the liquor flow increases 


so the lower velocity range is encroached and 
tends to disappear. 


2. Hold-up. 


Ve TRUE GAS VELOCITY ft/sec aT 15°C 


Fie. 9 Comparison of the various contributions to the 


hold-up. 
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The results show strong evidence for the 


Vain - Ye in COKE hypothesis that the hold-up of a packed system 
4 pepeeeaaayty waren Oe can be split into two portions, one of which is 
Q “aS flowing over the packing surface. 

FLOWING There remains one critical point to be 
| 0-04; An HOLD - UP PO emphasized, and that is the alteration in trend 
of Hp, n, K,, K, and K, as the coke size increases 
A | id until moe work can be done on larger coke 
4 WUE Y, VA YD sizes, but it is believed that the answer lies in the 

Uf GZ Y YY, fact that the water tends to stand out non-wettable 
° Za coke, and that, for the smallest size coke, the 
Ve in=3/ein COKE depth of the trickle running over the coke 
2 WATER RATE | WATER RATE | ; surface was of the same order as the size of the 
790 cuft/sqft hr 7976 cuft/sqft hr individual interstitial spaces. Thus it might 
3 A be expected, in retrospect, that the }’-}" coke 
ek size would show some anomalies, but it seems 
| Oost xan Ss unwise to say more until further evidence can be 

KW Lastly, to give a clear picture of the actual 
3 0-02 IFAT ys Ap > Hp value of the various contributions to the hold-up 
ML LA Hiquid flow and for one ten times 
oe greater. The figure requires little explanation. 
2 —— —— The equation represents the results satisfactorily 
© cuft/sqft hr 7-97 cuft,/sq ft he except for a minor discrepancy for the lowest 
F | liquid rates. It is noticed that the experimental 
points tend to follow a horizontal line for too 
long, and then to rise suddenly to fall on the 
a. theoretical curve at the higher values of the true 
| \ gas velocity. This discrepancy is attributed to 
siaiiadl inadequacy of liquor distribution which was 
$ . rectified when the gas velocity increased the 
3 —— > hold-up. It was not encountered at any of the 

7 to the Imperial Smelting Corporation Limi 
MLL VE) for permission to print this work. He also 
4 ° 4 


wishes to thank Mr. B. B. Baker and 
Mr. A. E. Tippinc for their aid in carrying 
out the experimental part, and Mr. P. A. T. 
Keeprine for advice in general. 
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NOTATION 


@ = cross-sectional area of flowing part of ky, kg = constants in pressure drop equations. 
the hold up (ft. ft.? cross-section of L = liquor rate (ft.3, hr ft.2 cross-section of 
packed column). packed column). 

b = perimeter of flow of flowing part of m,n, p, T = exponents. 
the hold-up (ft./ft.? cross-section of S = specific surface (total surface exposed 
packed column). to gas stream ft.2/ft.2 of packed 

d = characteristic dimentions of packing column). 
ane Vg =true gas velocity for air at 15°C, 

& = gravitational constant. based on free fractional cross-section 

H = total liquid hold-up (ft.*/ft.* of packed of column (fase). 

Hp = dynamic hold-up which will drain from 

the column (ft.5/ft.8 of packed column). 

ee pias mn). Subscript D for dry packing 

Hy = a contribution to hold-up at rest in the 
column due to liquid flow (ft.*/ft.* of .™ 

Hg = static hold-up which will not drain from p= demity. Subscripts 6, L and W 

refer to gas, liquid and water respectively 


the column (ft.4/ft.8 of packed column). 
A = pressure drop in gas stream (ft. of 7 = friction factor for a packed bed. 


irrigating liquor/ft. of column height). (f/k) = shape factor in equation (2). 
K,, Kg, Kg, K, = constants in hold-up equations. Re = Reynold’s number. 
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Optimum temperature gradients in tubular reactors—II* 
Numerical study 


O.ecu Bitoust and Neat R. Amundson 
University of Minnesota, Minneapolis 14, Minnesota 


( Received 28 May 1955) 


Abstract—In this paper extensive numerical calculations for the schematic reaction system 
A—» B-C are described in which the successive reactions may be of either first or second 


orders or both. From these computations the optimum temperature profiles and maximum 
yield of B for a given process time may be determined if complete data on the kinetics of the 
reaction are available. This work is a continuation of a previous paper. 


Résumé—Une étude numérique est présentée ci-dessous pour le systéme des réactions 
A—+ B-—->C ot les deux réactions consécutives peuvent étre indifféremment soit du premier, 
soit du second ordre. Ces calculs permettent de déterminer les profiles de température optimum 
et la production maximum du corps B dans un réacteur tubulaire si l'on dispose de la description 
cinétique compléte du systéme des deux réactions consécutives étudiées. Ce travail est la suite 


d'une publication antérieure. 


8. NUMERICAL StuDyY oF 
Two Consecutive Reactions 


Thus far the problem of determining the optimum 
temperature profile for a given reaction has been 
reduced to the solution of a system of simultaneous 
differential equations. This system is non-linear 
because of the non-linear dependence of the 
reaction rates on the temperature, and in many 
cases on concentrations. It would be impractical 
to search for formal solutions of such a system. 
Therefore the numerical solutions of the problem 
for several basic kinetic schemes will be presented 
below. These numerical solutions have been 
obtained by a simulation of the problem on an 
analogue computer (Reeves Electronic Analogue 
Computer, REAC) which is especially suitable for 
treatment of systems of non-linear differential 
equations. The results obtained on the computer 
are generally reproducible within 1%. They 
usually agree within 1% with results of exact 
solutions. The accuracy obtained depends upon 
the type of recording equipment used and on 
the accuracy of its calibration. The circuit 


* Part I—Chem. Eng. Sci., 1956 5 81-92. 


diagrams needed for the programming of the 
computer will not be included, but they can be 
supplied on request. 

To begin, a study of the shape of the tempera- 
ture profile was made for the case of two con- 
secutive first order reactions 


A+~+B-C 
with the use of equation (18) 
@T a—B;(dT\* dT 


Throughout this study the rate constants a (T) 
and b(T) have been assumed under the Arrhenius 
form 


Eiji 
a(T) = e® r) 


1 _ 1) 

b(T) =k, e* * 
where E, and E, are the activation energies of 
the two reactions and k,, and ky, are the initial 


values of the rate constants for the feed stream 
temperatures Ty. The notation 


+ Present address : Laboratorie Central des Poudres, Paris, France. 


¢ Notation as in Part 1, 
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= 


r= B/E, 
will be used. 
The temperature profile for a pair of consecutive 
reactions is always one of decreasing temperature 
as can be determined by an examination of 


dT na™ 
a = B(T) 


at least if the rate constants are of the usual 
exponential form. The influence of the ratio of 
activation energies, r, on the temperature profile 
has been investigated first, for the case 


n = 0-50, ty = 0-95, yo = 0-05. 


12010 
-|+ 30 —— 
' 


Je 25 


20 


TEMPERATURE 


00 10 20 
Time GROUP, 


Fic. 1. Study of the optimum temperature profile. Effect 
of activation energy ratio. 


A value y, = 0 was not taken because this leads 
to an infinite initial temperature gradient which 
is inconvenient on the computer. The results are 
presented in Fig. 1. It is seen that the temperature 
gradient in the first stage of the reaction becomes 
more severe when the ratio of activation energies, 
r, decreases. 

The second study involves the influence of the 
initial conditions, as expressed by the ratio of the 
initial rate constants, n, on the temperature 
profile as shown in Fig. 2. As might have been 


' 
Lad 


25 


TEMPERATURE GROUP 


(1) 
+ 060, «, = 095, y,* 00S 


10° 20 
time Group, k,t 


Fic. 2. Study of the optimum temperature profile. Effect 
of initial conditions. 


expected, the average initial temperature drop 
is higher for large values of n, that is to say for 
the case where initially the second reaction is 
faster than the first. It is to be expected in this 


Aw 
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Fic. 8. Study of the optimum temperature profile. Effect 
of feed stream composition. 
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Ye 


os 


o4 


02 


(1,1) 
x, #005, 


010 020 


case that more severe quenching is required. 
Finally, the influence of the initial amount of 
product B was examined by varying the ratio 
of 2/Yo with ry + Yo = 1. As shown in Fig. 3, 
the initial temperature drop was highest for low 
concentrations of product B. 

Up to now no account has been taken of the 
process time @ for which the calculated profile is 
the optimum one. Recalling that the zero 
function Z (t) has the form, equation (17), 


Z(t) =ax(t) — A(t) y(t) 
the value of Z(t) at time zero is 
Z (0) = — Yo. 
r 


Clearly Z (0) must be positive in all cases so that 
Fig. 3 must be used with some care. 


OPTIMUM YIELD 


Fic, 4. Optimum product yield against initial rate ratio. Both reactions first order. 


INITIAL RATE RATIO 


The general study of the case of two consecutive 
reactions, including the relation of initial con- 
ditions to process time, was conducted along 
somewhat different lines. First the feed con- 
centrations were chosen as 2) = 0-95, yg = 0-05, 
with no component C present. Although the 
detailed description of the setting of the problem 
on the computer will not be given, it is worth 
mentioning that a satisfactory circuit was obtained 
by taking the temperature as the computer time. 
The scheme adopted was A + B+ C, and the 
problem was treated for the following cases 


Order 
First reaction 


It was decided not to change the stoichiometric 
coefficients when both reactions are not of the 


Ye 


50/100 /o60 /040 /o10 | (1,2) 
110 090 070 050 030  anlideiinians 
00 L % 
010 020 050 100 200 


Fra. 5. Optimum product yield against initial rate ratio. First order reaction followed by a second order reaction. 
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Fic. 6. Optimum product yield against initial rate ratio, 


same order, because it was felt that no gain in 
generality would be obtained. 

The procedure adopted was to simulate the 
problem for a given value of the ratio of the 
activation energies, r, and for different values of 
the ratio of the initial rate constants, n, and to 
stop the solution when the zero function had 
crossed the value zero, thus indicating the 
optimum process time which corresponded to 
the initial conditions chosen. The results obtained 
in the four different arrangements of reaction 
orders were analyzed to give 
(1) The dependence of maximum yield of the 
product of interest, B, on the initial rate 
ratio, n (Figs. 4, 5, 6 and 7); 

The dependence of process time, @, on the 
initial rate ratio, n (Figs. 8, 9, 10 and 11); 


(2) 


OPTIMUM YIELD 


2.00 500 n 
INITIAL RATE RATIO 


Second order reaction followed by a first order reaction. 


100 


(3) The maximum temperature drop under 
optimum conditions as a function of the 
initial rate ratio, n (Figs. 12, 18, 14 and 15). 


These three studies were made for each of the 
four cases for different values of the ratio of 
activation energies, r. The results are presented 
in terms of dimensionless groups. Thus the 


temperature drop is represented by 7 


1 1 
T, 

where T, and 7, are the initial and final tem- 
peratures at time 0 and @, respectively. The 
reaction time is expressed by k,of and the initial 
conditions by n = ky/k,». For simplicity, the 
total number of reacting moles per unit volume is 
taken as one. This would necessitate a correction 
of the value of the second order constants in a 
case where this would not hold. The total 
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Fic. 7. Optimum product yield against initial rate ratio. Both reactions of second order. 
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Fic. 8. Optimum process time against initial rate ratio. Fic. 9. Optimum process time against initial rate ratio. 
Both reactions first order. First order reaction followed by a second order reaction. 


number of moles is assumed to remain constant. enumerated above, it seemed worthwhile to 

The effect of volume change may be included complement the theoretical methods of calcula- 

only in cases where it is.compatible with the tion by a thoroughly worked out case. It is felt 

kinetic representation assumed in this study. that the general trends which may be gathered 

Although the application of the charts to actua- from this study may at least lead to an apprecia- 

processes is thus limited by the restrictions tion of the problem in a given practical case. 
TiME GROUP TIME GROUP 
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Fie. 10. Optimum process time against initial rate ratio. Fic. 11. Optimum process time against initial rate ratio. 
Second order reaction followed by a first order reaction. Both reactions of second order. 
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Fic. 12. Overall temperature change against initial rate 
ratio. Both reactions of first order. 


Indeed, many practical processes involving con- 
secutive reactions may be schematized in such a 
way as to come under one of the cases here 
examined, so that a first study would show if it 
is economically advisable to try to increase 
reaction yield by controlling the temperature 
profile. 

The plots of maximum yield against initial rate 
ratio, n, show a general decrease with an increase 
in r and n. The maximum yield drops sharply 
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Fic. 14. Overall temperature change against initial rate 
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Fro. 18. Overall temperature change against initial rate 
ratio. First order reaction followed by a second order 
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to the initia? amount of B present in the feed 
mixture whenever the zero function at initial 
conditions is near zero. This occurrence, however, 
has little practical interest, since it is always 
advisable to have n as small as possible. It may 
also be noted that the peculiar aspect of the 
curves at the points where 
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ratio. Second order reaction followed by a first order Fic. 15. Overall temperature change against initial rate 


reaction, 


ratio. Both reactions of the second order. 
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is due mainly to our assumption that y, is not 
taken as zero. If y. = 0 the plots of yield vs. 
initial rate ratio, n, would have a regular sig- 
moida] shape. 

The trend with the different arrangements of 
reaction orders shows that high yields are mostly 
obtainable when the second reaction is of second 
order, which is easily explained, since in this case 
the destruction of the product of interest, B, 
would be smal] at the early stages of reaction 
because of the small amount of B present, and 
small during the later stages when the quenching 
is sufficient to cut down the reaction rate to a 
small fraction of its initial value. To visualize 
the effect of reaction order on the maximum 
yield, the maximum yield was plotted against 
the activation energy ratio, r, for the four 
different cases examined, for the value n = 1 of 
the initial rate ratio in Fig. 16. This figure 
actually represents cross-sections through n = 1 
of the plots of Figs. 4, 5, 6 and 7. 

A fact which may appear surprising at first 
glance is that the optimum temperature gradient 
is always negative whatever the value of the 
activation energy ratio, r. While the physical 
soundness of this fact is obvious for r < 1, since 
a temperature decrease in this case favours the 
formation of B, it is less apparent for r > 1. 
Indeed, when the activation energy of the second 
reaction is smaller than that of the first, a decrease 
in temperature is unfavourable to the ratio of 
the reaction rates, i.e. decreases the rate of 


ACTIVATION ENERGY RATIO 


Fie. 16. Maximum yield as a function of activation energy 
ratio for n = 1, for different arrangements of reaction 
orders. 


formation of B by a larger factor than the rate 
of destruction of B. Therefore, it might be 
expected that, comparatively speaking, a reaction 
with r > 1 would give a smaller yield under 
optimum temperature than a reaction with 
r <1. This conclusion must remain qualitative 
however, since such a comparison between pro- 
cesses of essentially different characteristics is 
rather difficult to establish on a sound basis. Such 
a qualitative trend is noticeable from the plots 
of optimum yield against initial rate ratio, and 
the comparison wil] be further developed in a 
numerical example. 

From the plots of optimum process time 
against the initial rate ratio, n, it appears that 
the value of the group k,.@ increases with 
decreasing n and increasing r. Further, it is 
always bigger when the second reaction is of 
second order. A conclusion with respect to 
process time, 8, instead of the dimensionless 
group k,. 6, is more difficult to achieve in all 
generality. An example is presented at the end 
of this chapter which helps to. clarify this 
relationship. The plots of overall temperature 
drop against initial rate ratio indicate that 
generally the temperature drop increases with 
the activation energy ratio, but is practically 
independent of it for small values of n. The 
temperature drop is more severe when the 
second reaction is of second order. 


9. APPLICATIONS AND 
NUMERICAL EXAMPLES 


The graphs presented thus far express the process 


time the overall temperature drop 
(7 7): and the maximum yield y,,, as a 
function of the initial rate ratio, n = kgo/kyo. The 
practical problem will usually occur as follows : 
determine the maximum yield, y,,, the overall 
temperature drop and the initial temperature, 
under optimum conditions, for a given process 
time, 

The rate constants, a (7) and 6 (7) are usually 
given in the form 


a(T) =k, (T) = 
b(T) = k,(T) = 
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The ratio n(7T) of these rate constants is 
+ (Bi — Bs 
n(T) = Pe 
1 
Eliminating the temperature, T, between k, (T) 
and n(T), one gets, 
k,(T) = py" 
Let 
a= 
then 


r 


k,(T) = 


The values of « and r may be calculated from 
the data of the problem, and the process time, 
@, is given in advance. Now 


(39) 


(40) 


may be plotted as a function of n on the general 
plot giving the process time group ky, @ against 
initial rate ratio n. On the log-log scale equation 
(40) plots as a straight line. The slopes of these 
lines depend upon r only, and are given together 
with the process time relationship for the case 
of two consecutive first order reactions on Fig. 17 
for r <1 and Fig. 18 for r > 1. 
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Fie. 17. Graphical determination of optimum initial con- 
ditions. Two consecutive first order reactions with r < 1. 
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Fic. 18. Graphical determination of optimum initial con- 

ditions. Two consecutive first order reactions with r > 1. 


As shown in Fig. 19, the line (1) represented 
by equation (40) intersects the process time 
group curve (2) at a point 4. The ordinate of A 
yields the value ky, of the rate constant of the 
first reaction at the initial temperature. The 
initial temperature is thus determined. The 
abscissa of the intersection 4 is the optimum 
initial rate ratio corresponding to the given value 
6 of the process time. This value of the initial 
rate ratio, when used on other plots will give the 
maximum yield, yyy, and the overall temperature 


change, corresponding to the 
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Graphical determination of optimum initial 
reaction rate ratio. 


Fie. 19. 
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given value @ of the process time. The problem 
is thus solved. Actually, an economic study will 
require the solution for a whole range of values 
of the process time @. Two numerical examples 
will now be considered. 

Example A—Case of two consecutive first order 
reactions where the activation energy ratio, r, is 
less than one, ie. Ey < By. 

The numerical values chosen are 


E, = 18,000 cal /mole 

E, = 30,000 cal /mole 

R = 2 cal /mol °K 

Pp, = 0-585 x 10" min 

P, = 0-461 x 10° 

6 = 10 min 

Ty = 0-95 mole litre 

Yo = 0-05 mole litre 

2% = 0 mole /litre 
Under these conditions 


k, = 0-585 x 7 
15 


,000 
k, = 0-461 x 10%¢@ 7 
and 


k 6,000 


— 0.868 x We 7 (41) 
ky 
k, = 21 


Construction of the straight line which represents 
this last relation on the graph of k,. @ vs. n given 
in Fig. 17 yields the initial condition 


n = 0-50 
This value, introduced into equation (41) yields 
the optimum initial temperature 

T, = 860°K 
From Fig. 4 which gives the maximum yield of 


component B for different values of the initial 
rate ratio, n, there results 


Ym, = 0-685 mole /litre 


From Fig. 18 with r = 0.60 the overall tempera- 
ture drop may be determined from 


or the drop in temperature is 24°K. Fig. 20 gives 
the concentration and temperature profile for this 
casel It is seen that most of the temperature 
change occurs in the first two-tenths of the 
process time. 
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Fro, 20. Optimum temperature and concentration profiles. 
Example A. 


MINUTES 


It is interesting to compare the yield under 
optimum conditions to the maximum isothermal 
yield. To obtain the latter, equation (42) below 
is solved for n,the group k, 8 being replaced by 
the expression 21 n**, There results 


n = 0-1852 
which corresponds to the temperature 
T = 340°K 


for the optimum isothermal process. The maxi- 
mum yield obtainable under isothermal con- 
ditions is 
Ym = 0-629 mole /litre 

Therefore in this case the optimum process 
gives a yield roughly 10% higher than the best 
isothermal process. If the reactions are exo- 
thermic it is expected that the adiabatic yield 
would be even less than the isothermal yield for 
the same initial conditions. It is interesting to 
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study the characteristics of the optimum opera- 
tion for different values of the process time. 
Such a study is presented in Fig. 21. The yield 
under optimum operation has been plotted against 
process time in a range extending from one to 
100 minutes. This yield may be compared with 
the optimum isothermal yield which is on the 


same graph. It is seen that the optimum yield 
increases with the process time. The initial and 
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ric. 21. Optimum yield and temperature requirements 


against process time. Example A. 


final temperatures and also the optimum isother- 
mal temperature have been plotted also. As the 
process time is increased the overall temperature 
drop tends to decrease. The optimum isothermal 
temperature lies between the extreme tempera- 
tures of the optimum process and closer to the 
final temperature. 

Example B—Case of two consecutive first order 
reactions where the activation energy ratio, r, is 
greater than one, i.e. E, > Ey. 

The numerical values chosen are 

E, = 30,000 cal /mole 
FE, = 20,000 cal/mole 
Pp, = 2-02 x 10 min" 
= 0-547 x 10" min” 
R = 2 cal /mole °K 

@ = 10 min 

yg = 0-95 mole /litre 

Yo = 0-05 mole litre 

= 0-00 mole /litre 
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The calculations are carried out essentially as 
in the preceding example. The results are 
Ym 0-465 mole litre 
T, = 370°K 
T, = 382.5°K 
Temp. drop = 37.5°K 
Best isothermal temp. = 340°K 
Best isothermal yield — 0-440 


Thus the gain in yield is very small indeed. A 
whole series of calculations has been carried out 
for different values of the process time, @, as in 
the preceding example, and the results are pre- 
sented in Fig. 22. It is seen that the yield in- 
creases when process time decreases, whereas the 


TEMPERATURE, 
a 


CONCENTRATION, 


to 20 50 20 30 
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Fic, 22. Optimum yield and temperature requirements 
against process time. Example B. 


inverse was true in the preceding case. The 
temperature drop and initial temperature also 
increase when the process time decreases. Thus 
when E, < E, the production of a high yield is 
limited by the reactor length, and when E, < E, 
by the heat transfer requirements. Similar con- 
clusions subsist when the two consecutive 
reactions are not first order but follow other mass 
action laws. 


10. IsoTHERMAL OPTIMA 


In section 2 the question of isothermal optima 
was examined for simple reactions. The purpose 
of this section is to compute the optimum 
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isothermal temperature for two first order con- 
secutive reactions 


Aw+~B-C 


the kinetic scheme being 


dt 
dy 


=> —k,« 


—kyy. 
The solution at constant temperature is 


k, 
—ket 
= Yoe + k, (e ). 
It is easily shown that the maximum of y occurs 
at a time 


Xo 
n n (Zp + Yo — NYo) 


and its value is 


™ [N(% + Yo — Yo) 


It is now desired to find the temperature at which 
¥» Will be a maximum for a fixed value of process 
time @. by 
dT 
and letting hp st k, be of Arrhenius form, the 
condition 
1 — NY 
x 
(1 — yo) (1 — 7) = eX 


r 


1+ 


(42) 


where X = (n —1)k,@ results. This trans- 
cendental equation defines n as a function of 6, 
that is to say, the optimum temperature as a 
function of the process time. Similar expressions 
are easy to devise for other kinetic schemes. 


ll. Orner PROBLEMS 
An examination of the zero function 


Z(t) =2(t) —A(y(t), A() = 2? 


a;(T)’ 


Equation (17), for the scheme A + B > C, each 
reaction being of the first order, shows that 
Z(@)=0 is the condition that the rate of 
formation of B at the outlet of the reactor (at 
t = 6) be a maximum with respect to tempera- 
ture. 

That is 


= way (T) 
=0 


— ybr(T) 


It has-frequently been suggested that such a 
condition holds at each point of the reactor. This 
is equivalent to the assumption that the yield of 
B will be a maximum at t = @ if the rate is at a 
temperature maximum at each point along the 
reactor. There is some foundation for this since 
for a simple reaction, section 2, the rate is chosen 
as a maximum at each point. However, if the 
condition 

by (T) E,k, 

y a7(T) E,k, 
in Arrhenius form, is substituted into the kinetic 


expression, equations (8 and 9), and equation (9) 
is divided by equation (8), there results 


The (E, — E,)/E, 
Thus : 


E, — E, 
_ = 
In example A this leads to 

y = 0.05 + (0.95 — 2) 


The maximum yield possible occurs when r=0 
and then y = 0-43 mole/litre. This yield is con- 
siderably smaller than either the optimum 
process or the optimum isothermal process. 
There are many other alternative schemes for 
the determination of optimum temperatures. 
These will not be exploited in full detail here but 


only suggested. 
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One might divide a reactor into two sections, 
the first section having a length fraction « of 
the total. The two sections might then be held 
at two constant temperatures 7, and T,, respect- 
ively. Thus the yield would be a function for 
given initial concentrations of three variables 


y =f (a, T,, T). 


The maximum yield y,, would then 
necessarily when 


dy oy 

and these three conditions could be used to 
determine a, 7,, and T,. These computations 
may be carried out using equation (11), for in 
this case the integrals contain integrands which are 
constants, different to be sure, over each section 
of the reactor. Three transcendental equations 
are obtained which must be solved simultaneously. 
This is not difficult to carry out, but it is tedious. 
In this way one can decide what temperatures 
the two isothermal sections should have, and also 
the relative reactor lengths for each section. It is 
clear that this may be generalized to any number 
of reactor sections in the same way. 

A variation on the above is the following. One 
might suppose that the sections of the reactor are 
not isothermal but, rather, adiabatic. The reactor 
could then be divided into two or more sections, 
each section having an intercooler. The optimum 
length of section and optimum intercooler tem- 
perature could then be determined in the same 
way as above. This is a more difficult problem, 
but can be carried out on modern computing 
machinery. 

Finally it is noticed that the optimum tem- 
perature profile appears to have the general shape 
of an exponential function, a parabola, or a 
higher degree parabola, and this appears to be 
quite generally true. If one supposes the profile 
may be expressed in the form 


occur 


=0 


T=a,+a,t+a,0@ 


with a,, a, and a, as constants, or as 
T=a,e" +a, 


then the yield is a function of a,, ag, a, and it may 
be maximized with respect to these constants. 
Such problems as this are almost as difficult to 
carry out mathematically as the determination 
of the optimum profile itself. 

From a practical point of view the impression 
of an optimum profile on a reactor might be 
difficult to carry out from the heat transfer point 
of view, and also the point control of temperatures 
might be a near impossibility. However, it seems 
clear that some scheme can almost always be 
devised for the approximation of optimum pro- 
files, and any mathematical study would have to 
be supplemented with more than a little engineer- 
ing expreience and art. It is only hoped that 
such a mathematical study may serve as a guide 
in a rational design, not as a law. 


CONCLUSIONS 


As it was pointed out in section 8, no rigorous 
proofs have been included in these papers in 
order to determine whether the quantities found 
are real maxima. These have been omitted. It 
is to be noted that the gradient of temperature 
must always be negative whether r > 1 or r < 1. 
In the latter case it is intuitively clear that the 
gradient must be a falling one. For r > 1 the 
situation is not quite so clear, and one’s intuition 
does not lead to such definite conclusions. The 
optimum isothermal yield is not much less than 
the optimum yield, and the difference is so little 
that small errors in computation either by hand 
or on the REAC might switch the results. In 
short, it is the feeling of the authors that more 
work on the case r > 1 is required. It is possible 
that the method employed in the paper does not 
give the physical solution when r > 1. 
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Die periodische Entstehung von Gasblasen an Dusen 
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Abstract—The volumes of yas bubbles, developed at small gas velocities on vertical, circular 
nozzles in liquids, were measured by using a stroboscope. Along with changing the gas velocity, 
the effects of variation of density, surface tension and viscosity of the liquid as well as the inner 
diameter and wall thickness of the nozzles were investigated. The formation of the gas bubbles 
in liquids of high and low viscosity was found to follow two different mechanisms. Two 
empirical equations based on the results obtained for the two ranges of viscosities were developed 
which permit calculation of the volumes of gas bubbles formed under these variable conditions. 


Résumé—A l'aide d'un stroboscope, les auteurs mesurent le volume des bulles gazeuses 
produites dans les liquides avec des vitesses gazeuses petites sur des ajutages verticaux ct 
5 circulaires. Tout en faisant varier la vitesse du gaz, ils étudient les effets du changement de densité, 
956 de tension superficiclle et de viscosité du liquide ainsi que leffet du diamétre intérieur et de 
lépaisseur des parois des ajutages. 

Ils trouvent que la formation de bulles gazeuses dans les liquides de grande et de faible 
viscosité suit deux mécanismes, Ils développent deux équations empiriques basées sur les 
résultats obtenus pour deux séries de viscosités, équations qui permettent de calculer les volumes 
des bulles gazeuses formées dans des conditions variables. 


INUALT EKINLEITUNG 


In einer kiirzlich erschienenen Arbeit wurde iber 
die Entstehung von Gasblasen an nach oben 
gerichteten, kreisférmigen Diisen berichtet [1]. Es 
wurde dabei der Mechanismus der Blasenbildung 
und -ablésung untersucht. Experimentelle Grund- 
lagen der Betrachtungen waren eigene Messungen 
am System Luft-Wasser sowie einige in der Litera- 
tur vorliegende Resultate, die sich auf andere 
Systeme bezogen. Solche Messungen, die den 
Einfluss der Flissigkeitseigenschaften auf die 
Blasenbildung zum Gegenstand haben — ein Ein- 
fluss des Gases ist im hier betrachteten Durchsatz- 
bereich nicht festzustellen — sind nur liickenhaft 
durchgefiihrt worden. In der vorliegenden Arbeit 
wird der Einfluss einer systematischen Variation 
verschiedener Eigenschaften der Flissigkeit auf 
die Blasengrésse im Bereich kleinerer Gasdurch- 
sitze, wo periodische Blasenbildung erfolgt, 
untersucht und eine Deutung der Ergebnisse 
gegeben, die sich an die friher mitgeteilten 
Vorstellungen anschliesst. Diese Vorstellungen 
sind hier, soweit es notwendig erschien, kurz 


skizziert. 


Leitet man durch eine nach oben gerichtete, 
kresiformige Diise (Gréssenordnung des Durch- 
messers : Millimeter) einen konstanten Gasstrom 
in eine Fliissigkeit, so bilden sich bei massigen 
Durschatzen (Gréssenordnung : einige cem/sec) 
streng periodisch einzelne Blasen, die nachein- 
ander aufsteigen. Steigerung des Durchsatzes hat 
zunehmende Unregelmassigkeiten in Blasenbildung 
und -aufstieg zur Folge, bis schliesslich bei hohen 
Durchsatzen (Gréssenordnung : 10 bis 100 cem /- 
sec) keine einzelnen definierten Blasen mehr an 
der Diisenmiindung entstehen, sondern ein 
Gasstrahl in die Fliissigkeit schiesst und dort in 
Blasen verschiedener Grésse zerfallt. Diese 
beiden Arten der Gaseinleitung kann man als 
“ Blasengasen ” und als “ Strahigasen ” bezeich- 
nen. In der vorliegenden Arbeit wird nur das 
Blasengasen, also der Bereich kleinerer Durch- 
sitze mit periodischer Blasenbildung an der 
Diisenmiindung behandelt. 

Damit die Messungen der Abhangigkcit der 
Blasengrésse von Diisenradius, Gasdurchsatz und 
F lissigkeitseigenschaften reproduzicrbar ausfallen 
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Apps 1. 


Apparatur zur Blasenerzeugung mit Pressiuft-zuleitung a, Ausgleichsgefiiss ¢ mit Entliftungshahn b, 


Stromungsmesser d mit auswechselbarer Kapillare ¢, Manometer f, Quetschhahn g, Diisenrohr h und Kivette i. 


und quantitativ gedeutet werden kénnen, missen 
einige Voraussetzungen erfillt sein. Einmal sollen 
die Blasen bei der Bildung auf der Innenkante der 
Diuse aufsitzen (vgl. die Abbn. 11 und 12), womit 


der Innenradius der Duse kennzeichnender 
Diusenradius wird (im Gegensatz etwa zu Abtropf- 
wo der Aussendurchmesser der Diise 
wesentlich ist). Weiter soll der Gasstrom durch 
den Dusenmindungsquerschnitt zeitlich konstant 
sein- unabhangig vom Bildungsstedium der 
Blase. Dann erweist sich bei héheren Viskositaten 
die Disenform als wichtig, die sich tber das 
Stromungsverhalten der Flissigkeit beim Ablé- 
sungsprozess auf die Blasengrésse auswirkt; sie 
muss also bei allen Diisen vergleichbar sein. Der 
praktisch wichtige Fall, dass die Gaszufuhr nicht 
konstant ist, sondern im Rhytmus der Blasen- 
bildung schwankt (etwa bei Siebplatten), lasst 
sich im Anschluss an die bei konstantem Gasstrom 
gewonnenen Ergebnisse erfassen, wie in der 
Arbeit 1) ausgefihrt ist. 


dusen, 


VERSUCHSANOKDNUNG 
DuRCHFUHRUNG DER 

Die Gasblasen wurden in der Apparatur der Abb 1 erzeugt. 
Aus der Pressluftleitung a wurde Luft in das 2 Liter grosse 
Ausgleichsgefiass c geleitet, wo durch den eng ausgezogenen 
Entliftungshahn b eine Feinregulierung des Luftstromes 
erfolgte. Cher den Differenzdruckstromungsmesser d mit 


UND 
MESSUNGEN 


auuswechselbarer Kapillare e und das Manometer f gelangte 
die Luft in das Disenrohr h, das in die Kiivette i von 6X8 
qem Grundfliche und 16 cm Hohe eintauchte, der Flissig- 
keitsspiegel befand sich im allgemeinen 4c¢m oberhalb 
der Diisenmiindung. Die Kiivette stand in einem Glas- 
behalter, durch den Wasser von 20°C aus einem Ultra- 
thermostaten geleitet wurde, sodass stets bei gleicher 
Flussigkeitstemperatur gemessen wurde. Der Quetschhahn 
g wurde so eingestellt, dass f mindestens 20 cm WS 
CUberdruck anzeigte. Damit war eine konstante Gaszufubr 
zur Diise gewihrieistet, da dann die stérenden Druck- 
schwankungen an der Disenmiindung bei der Blasenbildung 
klein gegen den Druckabfall am Quetschhahn waren. 

Die Diisen h bestanden aus doppelt geobogenen Glas- 
rohren verschiedenen Durchmessers, die mit einem Stick 
Schlauch ausserhalb der Fliissigkeit an die Luftzuleitung 
angeschlossen werden konnten. Die Diisenmiindungs- 


Asa. 2. Diise mit Blechring zur Erweiterung des Aussen- 
durchimessers. 
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flichen waren plangeschliffen, die Riander waren scharf- 
kantig. Der Aussendurchmesser sollte stets gross gegen 
den Innendurchmesser sein, um vergleichbare Stromungs- 
verhiltnisse bei der Blasenbildung und — abschniirung zu 
schaffen. Bei den Diisen mit kleinem Innendurchmessser, 
die aus Kapillarglas bestanden, war das ohne weiteres der 
Fall, bei den grésseren Diisen wurden nach Abb 2. Blech- 
scheiben um den Aussenrand gesetzt. Messungen ohne 
tliese Blechscheiben wurden ebenfalls durchgefiihrt, sie 
ergaben bei héheren Viskosilaten (etwa oberhalb 20 cPs) 
merklich andere Resultate. Tabelle 1 bringt die Innen- 
radien r, der benutzten Diisen sowie die Verhiiltnisse von 
Aussenradien r, zu Innenradien, 


Tabelle 1 


r, (in em) 
breiter Rand 
0,035 10,2 
0,059 5.9 
0,084 44 
O117 4) 2.8 
0,146 6,2 3.0 
0,195 60 3,1 


Das Aufsitzen der Blase auf der Diisenninenkante nach 
den Abbn. 11 und 12 wurde dadurch erreicht, dass die 
Diisenmiindung griindlich mit Chromschwefelsdure 
gereinigt wurde und somit benetzbar war, wahrend dic 
Innenwandung der Diise bis etwa 2cm unterhalb der 
Miindung leicht mit Siliconefett gefettet wurde, was cine 
Benetzung verhinderte. 

Die Messung des Blasenvolumesn J erfolgte strobos- 
kopisch iiber die Bestimmung der Bildungsfrequenz f. Da 
der Gasdurchsatz |" (Volumen /Zeit) am Stromungsmessecr 
abgelesen werden konnte, ergab sich }° nach 

vers (1) 
Als Stroboskop wurde eine von einem gut regelbaren 
Elektromotor getriebene Schlitzscheibe mit einem Schlitz 
benutzt. Sie war iber ein Zahnrad mit einem Schlagwerk 
gekoppelt, das bei 20 Umdrehungen einmal anschlug 
(Nachpriifung mit einer Glimmlampe). Zur Messung 
wurde die Kiivette von hinten mit einer 100 Watt-Lampe 
beleuchtet, der Motor so einreguliert, dass die Blasen bei 
Betrachtung durch die Schlitzscheibe stehen zu bleiben 
schienen und die Zeit fur 10-20 Glockenschlage abgestoppt. 
woraus sich leicht f ergab. Frequenzen unter 4 sec 
wurden durch direktes Abzahlen der enstehenden Blasen 
bestimmt. Der Fehler. mit dem das so bestimmte Blasen- 
volumen behaftet ist, betragt etwa + 5%. 

Zur Aufklarung des Blasenbildungsmechanismus wurden 
einige Photos gemacht. Dazu wurde eine Leica K mit 
Objektiv Focotar f = 5em 1: 4,5 und Wechselschlitten 
fiir Nahaufnahmen benutzt, die mit einem Elektronenblitz 
gekoppelt war. Der Blitzreflektor, der in Aufnahmerich- 
tung hinter der Kiivette stand, war mit einem weissen 
Blatt Papier abgeschirmt (Helifeldbeleuchtung). 


Die Bestimmung der Viskositat der untersuchten 
Fliissigkeit erfolgte in einem Ubbelohde-Viskosimeter mit 
hangendem Kugelniveau, die Messung der Oberflachen- 
spannung wurde nach der Blasendruckmethode durch- 
gefiihrt, die Dichte wurde ariiometrisch bestimmt. Alle 
Messungen wurden bei der Versuchstemperatur von 20°C 
durchgefiihrt. Von den genannten Fliissigkeitseigen- 
schaften wurde nach Moéglichkeit stets eine unter Kon- 
stanthaltung der anderen variiert. Es wurden dabei 
folgende Bereiche erfasst : 

Dichte p : 0.7 1.8 g/com 

Oberfliachenspannung o : — 73 dyn/em 

Viskositat » : 0,4 800 cPs 
Durchweg wurde fiir jede Fliissigkeit und jede der in 
Tabelle 1 aufgefiihrten Diisen das Blasenvolumen in 
Abhangigkeit vom Durchsatz bestimmt. 


Der VORGANG DER BLASENBILDUNG 


Im folgenden sind zusammenfassend die wichtig- 
sten Ergebnisse der Arbeit [1] referiert, dic 
Ausgangspunkt der vorliegenden Untersuchung 
waren ; von einer Wiedergabe der dort referierten 
Literaturstellen ist hier abgesehen worden. 


Die quasistatische Blasenbildung 


Drosselt man den Gasdurchsatz sehr stark, dann 
bildet sich die Blase so langsam, dass dynamische 
Effekte bei der Aufwélbung mit Sicherheit 
ausgeschlossen werden kénnen : die Blase hat in 
jedem Augenblick die nach der Kapillaritats- 
theorie zu erwartende Gestalt. Danach stehen in 
jedem Oberflachenelement der Blase der Gasinnen- 
druck, der hydrostatische Druck der Flissigkeit 
und der Kapillardruck, der von der Oberflachen- 
kriimmung herrihrt, miteinander im Gleichge- 
wicht. Aus dieser Gleichgewichtsbedingung 
ergibt sich cine Differentialgleichung fir die 
Gestalt von Blasenoberflachen, die approximativ 
losbar ist. Aus der Mannigfaltigkeit der Lésungen 
werden im hier betrachteten Fall durch die 
Randbedingung: Aufsitzen der Blase auf der 
Diiseninnenkante, d.h. bei festem Grundflachen- 
querschnitt, bestimmte Lésungen ausgewahlt, die 
nacheinander von der sich bildenden Blase 
realisiert werden. Schliesslich nimmt die (immer 
noch aufsitzende) Blase eine Gestalt an, die bei 
der gegebenen Grundflache ein maximales 
Volumen umschliesst. Lasst man weiter Gas 
zustrémen, sind keine mit der Kapillaritatsglei- 
chung vereinbaren Blasengestalten mehr méglich; 
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Ans. 3. Die maximalen Blasenbildungsvolumen fiir den 

quasistatischen Fall in Abhdngigkeit vom Disenradius 

(dimensionslos gemacht mit der Laplaceschen Konstante). 

exakter Verlauf, gestrichelt die lineare Naherung dafir; 
Verlauf fir ideale Kugelblasen. 


die Blase schnirt sich ab. Es zeigt sich, dass 
diese Abschniirung gerade so stattfindet, dass das 
Volumen der aufsteigenden Blase identisch ist 
mit dem theoretisch méglichen Maximalvolumen 
Vy, zwischen Disenmiindungsquerschnitt und 
Wandung der aufsitzenden Blase. WV, lasst sich 
nach der Kapillaritatstheorie als Funktion von 


r, und der Laplaceschen Konstanten a = Jf 


der Flissigkeit berechnen (o Oberflachenspannung 
Schwerebeschleunigung, p Flissigkeitsdichte). 
In Abb. 8 ist das dimensionslose Py, = Vy,/a* als 
Funktion von 7, = 1, ‘a (ebenfalls dimensionslos) 
aufgetragen. 

In der Literatur wird das quasistatische 
Blasenvolumen durchweg durch Gleichsetzen des 
Auftriebs einer als kugelférmig gedachten Blase 
mit der Haftkraft auf Grund der Oberflachen- 
spannung am Disenrand bestimmt. Neben der 
Kugelform der Blase wird bei diesem Ansatz 
implizit auch vorausgesetzt, dass der Diisenradius 
sehr klein gegen den Radius der sich abschniren- 
den Kugel ist. Beide Voraussetzungen sind nur 
in der Grenze fir verschwindende Disenradien 
erfullt, wie Abb. 8 zeigt, in der, wieder in dimen- 


sionsloser Form, das nach dem _ geschilderten 
Verfahren erhaltene Maximalvolumen als 
Funktion von 7, eingetragen ist ; diese Funktion 
ist linear. In guter Naherung kann fiir nicht zu 
grosse Diisen (7, < 0,7) auch My als lineare 
Funktion von 7, aufgefasst werden : 


2,5 7, (2) 


die Abweichungen vom wahren Verlauf betragen 
héchstens einige %. 


Die Messungen am System Luft-Wasser ergaben 
fiir eine bestimmte Diise eine lineare Abhangigkeit 
zwischen Blasenvolumen und Luftdurchsatz : 


V = Vy + (8) 


wo Wy, das der Diise (und der Flissigkeit) 
entsprechende quasistatische Blasenvolumen ist, 
und 7, eine von der Diise und den Flissigkeits- 
eigenschaften, nicht aber vom Durchsatz 
abhangige Grésse darstellt. Die Beziehung (3) 
lasst sich wie folgt deuten : 


Bei der Aufwélbung der Blase stellen sich auch 
hier die nach der Kapillaritatstheorie zu erwar- 
tenden Blasenformen ein, d.h. dynamische Effekte 
sind zu vernachlissigen. Ist Vy, erreicht, setzt 
der Abschniirvorgang ein. Da die Gaszufuhr 
konstant gehalten wird, strémt wahrend der 
Abschniirung weiter Gas in die Blase und zwar 
gerade das Volumen )’- T,, wenn T, die Zeit 
der Abschniirung bedeutet. Das Gesamtblasen- 
volumen V setzt sich somit additiv aus Vy und 
T,-V zusammen. Die Grisse T, in (8) ist also 
mit der Abschniirzeit der Blase identisch, und die 
Gleichung (8) sagt aus, dass die Abschniirzeit 
unabhangig vom Durchsatz ist. 7, betrigt einige 
hundertstel Sekunden. 

In der vorliegenden Arbeit soll nun untersucht 
werden, inwieweit dieser Blasenbildungsmecha- 
nismus fir andere Flissigkeiten Giltigkeit besitzt, 
wie, soweit das der Fall ist, JT, von den Flissig- 
keitseigenschaften und dem Disenradius abhangt, 
bzw. welche Abweichungen auftreten, und wie 
sie zu deuten sind. 
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Grenze des Blasengasens bei Durchsatzsteigerung 


Das Blasengasen wird bei hohen Durchsitzen 
vom Strahlgasen abgelést, der Ubergang zwischen 
beiden Arten der Blasenentstehung erstreckt sich 
uber einen weiten Durchsatzbereich. In diesem 
Bereich findet zwar unmittelbar an der Diisen- 
mundung noch eine periodische Abschnirung der 
Blasen statt, doch ist die Grésse der so entstehen- 
den Blasen nicht identisch mit der der aufstei- 
genden. Einmal tritt-vor allem bei engen 
Disen - die Erscheinung der Blasenverschmelzung 
auf: unmittelbar oberhalb der Diisenmiindung 
verschmelzen zwei Blaen zu einer. Weiter macht 
sich- vor allem bei den weiten Diisen — dic 
Instabilitat von mittelgrossen Blasen [2] bemerk- 
bar: die abgeschniirten Blasen zerfallen beim 
Aufstieg. und zwar um so leichter, je geringer die 
stabilisierende Oberflachenspannung der Flissig- 
keit ist. 

Hier wird nur der Fall relativ kleiner Durchsitze 
behandelt, in dem keine Komplikationen der 
gennanten Art auftraten. Aus genauen Beo- 
bachtungen der Verschmelzungs- und Zerfalls- 
vorgange beim System Luft-Wasser und den 
Ergebnissen der vorliegenden Arbeit bei anderen 
Systemen lasst sich eine empirische Forme! fiir 
den oberen Grenzdurchzatz V, dieses Bereichs 
angeben. Danach ist : 


v 


010 


i 


I’, 80 or,?, hichstens jedoch =, (4) 


wobei V, in cem/sec, o in dyn/em und r, in em 
anzugeben ist. 


Dir VERSUCHSERGEBNISSE UND 
IHRE DEUTUNG 
Aus den Messungen ergab sich, dass der oben fur 
das System Luft-Wasser geschilderte Mechanis- 
mus der Blasenbildung auch fiir andere Flissig- 
keiten Giltigkeit besitzt, sofern ihre Viskositaten 
nicht zu gross werden; als Grenze kann man 
etwa 20 cPs annehmen. Es sei hier zunachst auf 
die Vorginge in diesem Bereich eingegangen. 


BLASENBILDUNG IN FLUSSIGKEITEN 

NIEDRIGER VISKOSITAT (7 < 20 cPs) 
Einfluss der Dichte 
Mit verschieden konzentrierten Schwefelsduren 
wurden Dichten oberhalb 1 g/ceem bis zu 1,80 
g ccm hergestellt. Dabei wachst die Viskositat 
bis auf 25 cPs. Als Vergleichslésung wurden 
deshalb neben Wasser Glycerin-Wasser-Gemische 
entsprechender Viskositat gewahlt, deren Dichte 
nur wenig oberhalb 1 g/ccm liegt —die Unter- 
schiede in der Oberflachenspannnug sind in allen 
Fillen nicht sehr betrachtlich. Abb. 4 zeigt fir 
eine Diise die Blasenvolumen fiir Schwefelsiure 


2 
icPs 


WASSER 
circ 25 
H, SO, 25 


— 


i 


2 (cc m/sec) 


Awa. 4 Blasenvolumen V in Abhingigkeit vom Durch-satz p fiir Flissigkeiten verschiedener Dichte. 
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der Dichte 1,80 g/cem, Wasser und cin Glycerin- 
Wasser-Gemisch in Abhangigkeit vom Durchsatz, 
die Stoffkonstanten der drei Substanzen sind mit 
aufgefihrt. In allen Fallen liegen dic Werte 
gut auf Geraden. Die Abschniirzeiten (nach 
Gleichung (3) die Steigung der Geraden) stimmen 
praktisch tiberein. Die Schnittpunkte der 
Geraden mit der Ordinate, nach (3) die quasista- 
tischen Maximalvolumen, sind identisch mit den 
nach der Kapillaritatstheorie berechneten (Abb. 3) 
Gleiche Bilder ergeben sich auch fur die anderen 
Diisen, sowie fiir dic Messreihen mit Sauren 
anderer Dichte. 

Ein Einfluss der Dichte auf den Blasenbil- 
dungsmechanismus und die Abschnirzecit besteht 
somit nicht, die Dichte geht nur tiber die Lapla- 
cesche Konstante in das quasistatische Maximal- 
volumen Vy, und damit nach (3) in das Blasen- 
volumen V cin, 

Es sei aber angemerkt, dass bei sehr viel 
héheren Dichten (Quecksilber, Metallschmelzen) 
cin weiter gehender Einfluss der Dichte auftreten 


kénnte. 


Einfluss der Oberflichenspannung 


Tabelle 2 gibt cine Zusammenstellung der 
Flissigkeiten annahernd gleicher Viskositat und 


gleicher Dichte, jedoch verschiedener Ober- 


flachenspannungen, an denen der Einfluss dieser 
letzten Grisse auf die Blasenbildung untersucht 
wurde. Die Werte des Blasenvolumens V in 
Abhangigkeit vom Durchsatz V’ liegen bei allen 
Flissigkeiten und allen Diisen recht gut auf 
Geraden, wie Abb.5 am Beispeil der Essigsiure 
zeigt. Die leichten Schlangelungen der Messkurven 
um die repriisentierenden Geraden, in der Abb. 5 


Tabelle 2 


Oberfliichen- 
spannung 
in dyn ‘cm 

725 
42.0 
35.0 
28,5 


Viskositit 
in 
1,00 
154 
O95 
1,22 


Dichte 
in g/eem 
1,00 
1,22 
0,98 
1040 


Wasser 
Ameisensiiure 
Pyridin 
Essigsaure 


besonders an der Kurve fiir die Diise von 0,1945 em 
Radius erkennbar, scheinen allerdings keine 
Messfehler zu Im allgemeinen sind die 
Abweichungen linearen Verlauf jedoch 
klein. Die Ordinatenabschnitte stimmen wieder 
mit dem berechneten quasistatischen Maximal- 
volumen Vy, tiberein. 

Abb. 6 zeigt die Abschnirzeiten T , in Abhangig- 
keit vom Diisenradius fiir die untersuchten 
Flissigkeiten. In Anbetracht der Tatsache, dass 


sein. 
vom 


© O1945 cm 
o « 
40-084)" 
x 0593 


2 (cm/sec) 


Aun. 5. Blasenvolumen Vin Abhangigkeit vom Durehsatz fiir Essigsiiure an verschiedenen Diisen, 
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° or 

Ans. 6. Abschniirzeiten T , in Abhiingigkeit vom Disen- 

radius fir Flissigkeiten verschiedener Oberflachenspan- 
nung (Tabelle 2). 


die Bestimmung von 7’, als Anstieg mit relativ 
grosser Ungenauigkeit verbunden ist, liegen die 
Werte befriedigend in der Nahe eines Kurven- 
zuges, es ist dabei kein gesetzmassiger Einfluss 
der Oberflachenspannung auf die Lage zu dieser 
Kurve zu erkennen. Besser als durch eine 
friher verwendete lineare Beziehung lasst sich 
der Verlauf der Kurve durch eine Wurzelbeziehung 
wiedergeben : 
T = 0,052 4/7, (5) 
wobei 7, in sec, r, in cm einzusetzen ist. 
Durch Zusatz oberflachenaktiver Substanzen zu 
Wasser lisst sich dessen Oberflachenspannung 
stark andern, ohne dass Dichte und Viskositat an 


(ccm) 


dieser Anderung teilnehmen. Es wurden deshalb 
Untersuchungen an einer 1% -igen wassrigen 
Lésung des im Handel erhaltlichen Mittels “ Pril ” 
vorgenommen, deren Oberflachenspannung 39 
dyn/em betrug. Die Kurven V in Abhangigkeit 
von V’ steigen dabei fiir kleine Durchstatze vom 
quasistatischen Wert ausgehend annahernd 
linear an, wobei die Steigung (Abschniirzeit) 
jedoch merklich grésser ist als bei Wasser. Bei 
héheren Durchsatzen knicken die Kurven ab und 
nihern sich asymptotisch der entsprechenden 
Geraden fiir Wasser. Dies Verhalten wird darauf 
zurickzufiihren sein, dass zur Einstellung der im 
statischen Versuch gemessenen Oberflachen- 
spannung eine gewisse Relaxationszeit erforder- 
lich ist, die bei grésseren Durchsatzen wegen der 
Schnelligkeit der Oberflachenneubildung nicht 
mehr zur Verfiigung steht. Méglicherweise lasst 
sich auf diese Erscheinung ein Verfahren zur 
zumindest rohen Abschatzung solcher Relaxa- 
tionszeiten aufbauen. Fir die vorliegende Unter- 
suchung waren die Messungen nicht brauchbar. 

Auf Grund der Abb. 6 lasst sich feststellen, dass 
auch die Oberflachenspannung ausser ihrer 
Wirksamkeit fiir das quasistatische Maximal- 
volumen keinen Einfluss auf Blasenbildung und 
Abschniirzeit hat. Es muss dabei vorausgesetzt 
werden, dass sich die Oberflachenspannung bei 
verschieden schnell erfolgender Blasenbildung 
nicht merklich andert. 


© 0-945 cm 
+ O70 
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Ass. 7. Blasenvolumen V in Abbangigkeit vom Durchsatz / fir m-Kresol an verschiedenen Diisen, 
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Einfluss der Viskositat 

Es wurden Fliissigkeiten annihernd  gleicher 
Dichte und Oberflaichenspannung, jedoch ver- 
schiedener Viskositat, wie sie in Tabelle 3 auf- 
gefiihrt sind, untersucht. Wieder lassen sich die 
Blasenvolumen in Abhangigkeit vom Gasdurch- 
satz durch eine lineare Funktion erfassen, wie 
Abb. 7 am Beispiel des m-Kresols zeigt. Abb. 8 


x WASSER 
@ m-KRESOL 


CS; 
© PYRIDIN 


iL 


Oo" 


Ass. 8. Abschniirzeiten T 4 in Abhingigkeit vom Diisen- 
radius fiir Fliissigkeiten verschiedener Viskositit (Tabelle 
3). 


Tabelle 3 


Obertliichen- 
spannung 
in dyn/em 
33,8 
33,0 
36,0 


Viskositat 
in 


Dichte 
in g/eem 
1,263 0,366 
0,98 
1,054 21 


Schwefelkohlenstoff 
Pyridin 
m-Kresol 


bringt die Abschnirzeiten in Abhaingigkeit vom 
Diisenradius, Der friher vermutete starke Einfluss 


der Viskositat auf dic Abschniirbewegung, die 
doch eine Scherbewegung darstellt, zeigt sich 
merkwiirdigerweise nicht, die Abschniirzeiten fiir 
Kresol liegen nur unwesentlich oberhalb der 
anderen Werte. Die Bezichung (5) hat, wie dic 
Abb. 8 zeigt, auch fiir die Abschniirzeiten bei 
Flissigkeiten verschiedener Viskositat Giltigkeit, 
solange diese Viskositat 20 cPs nicht merklich 
tiberschreitct. 

Viskositaten unter 20 haben also praktisch 
keinen Eintluss auf die Blasenbildung. 


Formelmdssige Darstellung der Ergebnisse 

Aus den Gleichungen (2), (3) und (5) ergibt 
sich als explizite Formel fiir das Volumen V (ccm) 
einer Gasblase, die in einer Fliissigkeit der 
Obertlichenspannung (dyn/cm), der Dichte p 
(g/eem) und der Viskositét < 20 cPs an einer 
Diise des Innenradius 7, (em) beim Gasdurchsatz 
(com /sec) entsteht : 


(6) 


V = 5,110" + 5,2- 10° 74/7, 


Voraugsesetzt wird dabei konstanter Gasstrom 
durch die Diise, Aufsitzen der entstchenden 
Blasen auf der Diiseninnenkante und cin Durch- 
satz <V,, wo, = 80072, hichstens jedoch 


5 Ein Einfluss der Form des Diisenrands (breit 


oder schmal) war zumindest aus Messungen an 
den in der Tabelle 1 aufgefiihrten Diisen nicht 
zu erkennen. 

In Tabelle 4 sind einige Beispiele fir gemessene 
und nach Gleichung (6) berechnete Blasenvolumen 


zusammengestellt. Es zeigt sich eine gute 


Tabelle 4 


dyn/em 
Wasser 72,5 


Ameisensiure 42,0 
33,8 
35.0 
28.5 
36.0 


Schwefelkohlenstoff 
Pyridin 

Fssigsiure 
m-Kresol 


Schwefelsiiure 58,0 


V exp. 
ecm 
0,0538 
0,1215 
0,0328 
0,0576 
0.0131 
0,0295 

0,014 
0,072 
0,0884 
0,0761 


V ver. 

com 
0,054 
0,122 
0.0316 
0,0134 
0,0279 
0,0146 
0,076 
O,0757 
0,0696 


(oe) 
a 
oo 
: || VOL. 
5 
1956 
° 
r, v 
g/ecm em ccm /sec 
1,00 0,084 1,0 
4,0 
1,22 0,059 1,06 
0,146 1,24 
1,263 0,035 0,58 
0,98 0,084 0,55 
1,049 0,059 0,318 
0,084 0,117 2,0 
0,195 1,37 
1,80 0059 
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Ans. 9. Blasenvolumen V in Abhangigkeit vom Durchsatz j’ fiir ein 
Glycerin-Wasser-Gemisch von 75 cPs an verschiedenen Diisen. 
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Ans. 10. Blasenvolumen in Abhingigkeit vom Durch-satz fiir 
ein Glycerin-Wasser-Gemisch von 600 cPs an verschiedenen Diisen. 
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Obereinstimmung ; selbst beim Kresol und der 
Schwefelsiure, deren Viskositaten schon oberhalb 
der 20 cPs-Grenze liegen, weichen gemessene und 
berechnete Werte im Mittel nur um 10% 
voneinander ab. 


BLASENBILDUNG IN FLUSSIGKEITEN 

ViskositTAt (» > 20cPs) 
Bei Flissigkeiten niedriger Viskositat crwiesen 
sich Oberflachenspannung und Dichte als nicht 
wesentlich fir den Abschniirprozess. Es kann 
angenommen werden, dass dies bei Flussigkeiten 
hoherer Viskositat mit dem tberragenden Einfluss 
der Viskositat auf den Blasenbildungsvorgang 
erst recht der Fall ist ; deshalb wurden in diesem 
Bercich Dichte und Oberflachenspannung nicht 
systematisch variiert. Bei hohen Viskositaten ist 
es auch schwer, Flussigkeiten sehr verschiedener 
Oberflachenspannungen zu finden. 

Die Ergebnisse zeigten eine Abhangigkeit von 
der Dusenform, d.h. von der Breite des Dusen- 
randes. Zunachst seien die Messungen an 
breitrandigen Diisen (Tab. 1) mitgeteilt und 
besprochen. 


Mesungen an breitrandigen Diisen und thre 
Deutung 


Als Flussigkeiten dienten 60-70% —ige Rohr- 
zuckerlésungen, die Viskositaten bis zu 230 cPs 
aufweisen, und Glycerin-Wasser-Gemische mit 
Zahigkeiten bis zu 700 ¢Ps. Die Dichten betrugen 
bei @n Zuckerlésungen 1,25-1,35 g cem, bei den 
Glycerin-Wasser-Mischungen 1,20-1,25 g cem, die 
Obertlachenspannungen 72-76 dyn bzw. 
64-66 dyn cm. 

Dic Abb. 9 zeigt fiir verschiedene Dusen des 
Blasenvolumen in Abhangigkcit vom Durchsatz 
fir cin Glycerin-Wasser-Gemisch von 75 cPs, dic 
Abb. 10 des Gleiche fir ein solches von 600 cPs, 
Man erkennt in allen Fallen ein starkes Ansteigen 
des Volumens mit dem Durchsatz fur sehr kleine 
Durchsatze, wobei fiir J’) +0 die berechneten 
quasistatischen Werte erreicht werden. Dann 


knicken bei etwa J’ = 0,2 com sec die Kurven 
ab und verlaufen als Geraden mit Steigungen, clic 
den bei Flissigkeiten mit Viskositaten unter 
20 «Ps auftretenden vergleichbar sind. Es gilt 
in diesem Bereich also : 
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(7) 
wo V,. der Schnittpunkt der Geraden mit der 
Ordinate, nicht mehr mit dem quasistatischen 
Blasenvolumen identisch ist. Eine Deutung 
dieser Erscheinung lasst sich auf Grund folgender 
Vorstellung geben : 


Die sofortige Ausbildung des quasistatischen 
Blasenvolumens auch bei hohen Durchsatzen, 
d.h. schnellen Flissigkeitsbewegungen oberhalb 
der Diise, wie sie bei den Viskositaten unter 
20 cPs erfolgt. tritt nur ein, wenn dynamische 
Effekte infolge der  Flissigkeitsbewegungen 
vernachlassigbar sind. Das ist nun bei den 
hohen Viskositaten nicht mehr der Fall. Steigt 
der Durchsatz vom Wert Null aus an. so bewirkt 
die innere Reibung der Flussigkeit bei der 
Aufwoélbung der Blase Deformationen gegenuber 
der quasistatischen Blasengestalt. Das Volumen, 
nach dessen Erreichen die Abschnurung einsetzt, 
wird dabei grésser als das quasistatische Volumen. 
Schliesslich wird bei weiterer Steigerung des 
Durchsatzes eine Art der Aufwolbung erreicht, 
bei der ein Minimum an Scherbewegungen auf- 
tritt (weitgehend Kugelform der sich aufwolben- 
den Blase; eventuell durch Sogwirkung der 
aufsteigenden Blasenvorgangerin nach oben 
verlangerte Formen). Das Volumen, nach dessen 
Erreichen hierbei die Abschnurung beginnt, ist 
V,. Diese Art der Blasenbildung findet sich dann 
auch bei allen hoheren Durchsatzen. Der lineare 
Verlauf der Abhangigkeit des Blasenvolumens 
vom Durchsatz erklart sich wieder wie bei kleinen 
Viskositaten dadurch, dass wahrend der Ab- 
schnirung, die unabhangig vom Durchsatz in der 
konstanten Zeit T, erfolgt, das Gasvolumen 
T,-V in die Blase einstromt. das sich zu Vy 
addiert. Im folgenden sei nur der Durchsatz 
besprochen, in dem Gleichung (7) gilt ; der Anlauf 
bereich ist schwer zu erfassen und wegen seincr 
kleinen Durchsatze (unter 0,2 ccm sec) praktisch 
bedeutungslos. 

In den Abbn. 11 und 12 sind zum Vergleich 
cinzelne Bildungsstadien von Blasen an der 
gleichen Duse und beim gleichen Durchsatz einmal 
in Wasser und zum anderen in einem Wasser- 
Glycerin-Gemisch von 600 cPs zusammengestellt. 
Sie illustrieren die obigen Ausfuhrungen. 


Ass. 11. Bildungsstadien einer Luftblase an einer Diise 
von 0,084 em Radius in Wasser. Durchsatz 0,65 cem/sec ; 
Blasenvolumen 0,046 ccm. 
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Ans. 12. Bildungsstadien ciner Luftblase an einer Diise 

von 0,084em Radius Glycerin-Wasser-Gemisch von 

Goo Durchsatz 0.65 cem sec; Blasenvolumen 
0,155 cem. 
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Die Werte fir V, aus den verschiedenen Mess- 
reihen sind in Abb. 13 in Abhangigkeit von der 


Ih 


Viskositat mit dem Diisenradius als Parameter ( 


aufgetragen. Wie die T,—Werte sind auch die 
V, ~— Werte nur relativ ungenau zu bestimmen, sie 
ordnen sich fiir die einzelnen Diisen aber befrie- 
digend um Geraden an; die Schnittpunkte der 
Geraden mit der Ordinate ergeben dann das 
quasistatische Blasenvolumen. Die Priifung der 


Abhangigkeit der Steigungen der Geraden vom 
Diusenradius ergibt cin Wurzelgesetz, sodass man 
zu der empirischen Formel gelangt : 


Vo = + 8,16-10% 4/r, 


(8) 


wo Vy und Vy in cem, » in cPs und r, in em 
einzusetzen sind. 

Zur Bestimmung der Abhingigkeit der 
Abschniirzeit T., vom Diisenradius r, und der 
Viskositat » sind in Abb. 14 die T., — Werte gegen 
r, als Argument mit » als Parameter aufgetragen. 
Fir gleiche Viskositaten liegen die Punkte 
recht gut auf Geraden. Nach einer schon friher 
mitgeteilten Vorstellung lasst sich das deuten, 
indem man 7’, in zwei Teilzeiten zerlegt : In der 
ersten, vom Radius unabhangigen, wird nach 


+ cm 
cm 
a cm 
© 0-084) cm 
xX cm 


Asp. 18. Blasenvolumen I'9, nach dessen Erreichung die 
Abschnérung einsetzt, in Abhiingigkeit von der Viskositat 
fir verschiedene Disen. 


° 


Ans. 14. Abschniirzeiten T , in Abhangigkeit vom Disen- 
radius fiir verschiedene Viskositaten. 


Erreichen von V, die Abschniirung eingeleitet, in 
der zweiten zieht sich der abschniirende Flissig- 
keitswulst mit annahernd konstanter Radial- 
geschwindigkeit zusammen, diese Teilzeit kann 
als proportional dem Radius angenommen werden. 
Fir die Radialgeschwindigkeit ergibt sich in 
Ubereinstimmung mit friheren, sich auf das 
System Luft-Wasser beziehenden Uberlegungen 
8-10 em sec, hier macht sich also kein wesent- 
licher Einfluss der Viskositat bemerkbar. 

Mit einer empirischen Formel lasst sich T, 
(sec) in Abhangigkeit von » (cPs) und r, (cm) 
darstellen durch : 


T , = 0,015 + 10° »'* + 0,097, (9) 


Einfluss der Diisenform 

Messungen an diinnwandigen Diisen zeigen 
Abweichungen gegeniiber Bestimmungen des 
Blasenvolumens an dickwandigen Dusen gleichen 
Innendurchmessers (Tab. 1). Abb. 15 zeigt das 
fiir die Diise von r, = 0,146 cm. Wahrend V, 
ibereinstimmt, unterscheiden sich die Steigungen 
der Geraden, d.h. die Abschniirzeiten T ,: bei den 
diinnwandigen Diisen ist 7, kleiner. Abb. 16 
bringt T, in Abhangigkeit vom Radius r, fir 
zwei verschiedene Viskositaten : mit wachsender 
Viskositat werden die Abweichungen zwischen 
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25 
Ass. 15. Blasenvolumen in Abhangigkeit vom Durch- 


satz J fir Disen gleichen Innendurchmessers aber ver- 
schiedenen Aussendurchmessers bei hohen Viskositéten. 


den Werten fiir dinn-—und dickwandige 
Diisen grésser. 

Die Erklarung fiir diese Erscheinung ergibt 
sich aus einer Betrachtung der Strémungsvor- 
ginge bei der Ablésung. Wenn die Blase auf- 
steigt und sich abschniirt, muss in den Raum 
oberhalb der Diise Flissigkeit einstromen. Kann 


das bei den dinnwandigen Diisen nach Abb. 17a 
der Fall ist, geht die Abschniirung rascher vor sich, 
als es bei den dickwandigen Diisen der Fall, ist wo 
die Flussigkeit nur von der Seite her nachstrémt. 
(Abb. 17b). 

Um diese Erscheinung zu erfassen, muss 
systematisch das Verhaltnis von Disennaussen- 


die Flissigkeit auch von unten nachfliessen, wie radius zu Disenninenradius von 1 bis zu 


l 
0-20 


ac) 


ous 


Ass. 16. Abschnirseiten T , in AbhAngigkeit vom Radjus 

fir Disen gleichen Innendurchmessers aber verschiedenen 

Aussendurchmessers (dick- und dinnwandige Diisen) bei 
twei Viskosit&ten. 


v 
(OIC Kw.) 
5 
200 ¢ Ps (DICKW) 
(Ummm) 
188 
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Die periodische Entstehnug von Gasblusen an Diisen 


variiert werden. Die oben mitgeteilten und 
diskutierten Werte entsprechen wohl annahernd 
dem Verhaltnis o. 


ja) (b) 


Ass, 17a und 17b. Flissigkeitsstrémung bei der Abschnii- 
rung an diinn- und dickwandigen Diisen. 


Formelmassige Darstellung der Ergebnisse 


Fiir dic Abhangigkeit des Blesenvolumens von 
den hier wichtigen Einflussgréssen wurde von 
A. und E. Pereruans [3] cine Formel 
vorgeschlagen, die nur als Grenzgesetz fiir sehr 
kleine Durchsatze Geltung beanspruchen kann, 


da sie den Durchsatz nicht enthalt. Aus den 
Bezichungen (7), (8) und (9) ergibt sich die 
allgemeinere Forme! : 


51-1072 + 3.16.10 + 
p 


(0.015 + 10% 0.09 (10) 


wobei die Oberflachenspannung o in dyn em, dic 
Dichte p in g cem, der Diisenradius r, in em. dic 
Viskositat » in cPs und der Durchsatz Vin 
cem sec einzusetzen sind. Die Formel gilt fir 
Durchsatze oberhalb 0,2-0.3 cem see und unter- 


halb 30 héchstens jedoch 5 ccm/sec, sowie 


fiir sehr breitrandige Dusen. Die Blasen miissen 
auf der Diseninnenkante aufsitzen, die Flissig- 
keit darf die Diisenwand nicht benetzen. 

In Tabelle 5 sind einige gemessene und nach (10) 
berechnete Blasenvolumen zusammengestellt, die 
Ubereinstimmung ist recht gut. 

Dem Leiter des Institutes, Herrn Prof. Dr. 
dinken wir fir wohlwollende Férderung 
der Arbeit und sein standiges Interessc. 


Tabelle 5 


dyn/em 
65 


cm 4 com 
0.0593 0,128 
0,286 
0,0593 0,174 
0,366 

0,0593 0.0588 
0,126 


ZUSAMMENSTELLUNG DER SYMBOLE 


- bedeutet, dass die Lingenpotenz 2 durch Division 
durch die entsprechende Potenz von a dimensionslos 
gemacht wurde 


a = Laplace’sche Konstante = Jz (em) 


Jf = Frequenz (1, sec) 

g Schwerebeschleunigung sec?) 
~ Aussenradius der Diise (em) 

= Innenradius der Diise (cm) 


‘a 


p Flissigkeitsdichte (g 
o = Oberflachenspannung (dyn/cm) 
Abschniirzeit (sec) 
= Blasenvolumen (cm*) 
J’ = Durehsatz (em? /sec) 
= scheinbares quasistatisches Blasenvolumen (cm?) 
= oberer Grenzdurchsatz 
= theoretisch mégliches Maximalvolumen (em?) 


= Viskositaét (cT’s) 
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5 
V ber. 
cPs com 

Glyc.-Wasser 430 1,25 0,122 

0,267 
610 64,5 1,25 0,170 

0,385 

Zuckerlos 130 75 1,32 0,057 
0,116 


J. M. Covursow and J. F. Ricnarpson: Chemical 
Engineering, Vol. II, 1055. Pergamon Press Ltd., London, 
and McGraw-Hill, New York. XVI + 588 pp. 60s. ($9). 
Tue subtitle of this volume is “ Unit Operations.” As 
a logical successor to Volume I, which was devoted to the 
basic concepts of fluid flow, heat and mass transfer, it 
contains the application of these concepts to the physical 
operations of material handling, mixing and separation. 
Naturally, the separation methods fill the greater part of 
this substantial book. 

The authors, who have themselves contributed valuable 
and original material to this field in the course of many 
years, and who both have a wide teaching experience, 
have subdivided the book into five main sections. As far 
as possible a systematic division has been attempted, 
mainly on the basis of the mechanism involved. So we 
have the following sections : 


E — Flow of fluids past particles (fixed beds and packed 

columns ; filtration ; the centrifuge), 

Relative motion between fluid and particles (particles, 

bubbles, and drops; sedimentation, fluidization, 

and conveying ; gas cleaning), 

G - Applications of mass transfer (leaching ; distillation ; 
gas absorption ; extraction), 

H - Evaporation, crystallization, and drying, 

I — Size reduction and classification of solids, mixing. 

The chapter on distillation is the longest one and contains 

nearly 100 pages. 

It is hardly possible to present a detailed discussion 
of the various chapters. The amount of information 
collected in this book is overwhelming, and one wonders 
whether im the future it will still be possible to have all 
unit operations together in one volume of reasonable size 
without the omission of a number of essential aspects. 
Moreover, different groups of operation may require a 
rather different treatment. Thus, according to whether 
the topic is already a “ science" or still an “ art,” one 
might analyse mechanisms and rate-determining factors 
or be merely descriptive. In between we may have the 
inevitable empirical correlations, multi-stage calculations, 
trial-and-error methods, etc., which in some textbooks 
tend to draw a curtain between the reader and the subject 
matter. 

In this book the authors have, as much as possible, 
stressed the mechanisms involved in the various operations, 
thereby frequently referring to the matcrial of Volume I. 
This starting point is well expressed by the final sentence 
of the authors’ introduction: “It is only by under- 
standing the mechanisms, and then the method of con- 
structing the equipment, that it will be possible to obtain 
a complete picture, including the cost of operation.” 
Of course, the fundamental aspects had to be amplified 
to a great extent by typical chemical-engineering calcula- 
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tions (multi-stage processes in Section E) and sometimes 
by much empiricism (e.g. in Section I). But it is exactly 
the authors’ point of view cited above which makes the 
book as a whole very valuable, not only to chemical- 
engineering students, but also to chemical engineers in 
practice. It is easy to read, thanks to a clear presentation 
and style. 

Since the book has only two authors, the various topics 
have been treated in a “ homogenous "' manner. For the 
reader, this is a great advantage which outweighs the 
advantages of complete coverage of all the material 
by a group of experts. Many literature referenoes conclude 
each chapter, including a number of very recent ones. 
The 121 problems at the end of the book will prove to 
be very useful. 

H. Kramers 


A. J. Younc: An Introduction to Process Control 


System Design, 1955. Longmans Green and Co. Ltd., 
London. XVII + 875 pp. 42s. 


A book by an acknowledged authority is always to 
be welcomed. This book is, as would be expected, 
competent and thorough. Its aim is that of “ interesting 
the reader who is not necessarily an instrument specialist 
but who is nevertheless concerned with automatic control 
as a tool or technique for increasing plant efficiency.” 

The presentation is very satisfactory, and the text is 
liberally interspersed with diagrams. The latter, though 
clear and adequate, are sometimes rather small. Minor 
errors only were noted, though exception must be taken 
to the use of “ calorific value of steam (p. 27) when 
enthalpy is meant. 

An introductory chapter on economic aspects gives 
powerful arguments for increased application of automatic 
control on economic grounds. A chapter on the closed 
loop is followed by five chapters devoted to the characteris- 
tics, controllability, and response of process plant, and 
six on controllers in theory and practice including descrip- 
tions of the essential mechanism of a wide range of 
industrial makes of controller. There follow chapters on 
measuring and transfer lags, transfer stages, valve charac- 
teristics, and aspects of control system design. Six 
appendices deal with topics such as control terminology, 
contro'ler adjustments, and the automatic factory. 

Discussing report writing, Kinxsatpe in one of his 
books quotes the advice given to a young lawyer, “ Tell 
the judge what you are going to tell him. Then tell him. 
Finally tell him what you told him.” This admirable 
method is followed here, the ideas of each chapter being 
nicely crystallized before the next is tackled. The book 
is thus extremely readable, and should fulfil the aim 


quoted above. 
D. M. Wrison 


